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sunriARY 


The chemical behavior of jet aircraft emissions In the atmosphere was 
studied using computer simulations of static well mixed air parcels. Emissions 
data contained In the literature Indicate that 85 to 95 percent of the hydro- 
carbons emitted from jet aircraft In the vicinity of airports are associated 
with the taxi-idle mode. The types and unounts of compounds present are 
listed below: 


Jm 


Percent Average 

(by volume) Carbon Number 


Paraffins 50-60 8 

Olefins 15-30 3 

Aromatics 7-30 8 

Aldehydes 1-17 2 


Nitrogen oxide emissions are associated primarily (80 to 90 percent) with the 
takeoff and cllmbout modes, and we estimate a 9:1 mole ratio of nitric oxide 
to nitrogen dioxide. Most significantly, the literature data show that the 
ratio of hydrocarbons to nitrogen oxides averages nearly 42 moles of carbon 
(C) per nxjle of nitrogen oxides (NO^). Observations at airports confirm this 
emissions ratio; those data average about 45 moles C per mole NO^. 


The emissions data were used as initial conditions for a series of 
computer simulations of photochemical smog formation in static air. The 
chemical kinetics mechanism used in these simulations was an updated version 
of the Hecht, Seinfeld, and Dodge (1974) mechanism. 


This kinetics mechanism contains certain parameters which are designed to 
account for hydrocarbon reactivity. These parameters were varied to simulate 
the reaction rate constants and average carbon numbers associated with the 
jet emissions. The roles of surface effects, variable light sources, NO/NO 2 







V 

ratio, continuous emissions, and untested mechanistic parameters were also 
assessed. 

The results of these calculations indicate that the present jet emissions 
are capable of producing oxidant by themselves. The hydrocarbon/NO^^ ratio of 
present jet aircraft emissions is much higher than that of automobiles. These 
two ratios appear to bracket the hydrocarbon/NO ratio that maximizes ozone 

A 

production. Hence an enhanced effect is seen in the simulations when jet 
exhaust emissions are mixed with automobile exhaust emissions. 
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I INTRODUCTION 
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Tlw fonnatlon of photochemical oxidant has long been known to be related 
to Interactions of two classes of prlnary pollutants* the hydrocarbons and 
the nitrogen oxides. These pollutants are normally associated with the use 
of fossil fuels: hydrocarbons from Incon^lete c(xitust1on or fuel leakage 
and nitrogen oxides from exposure of air to high temperatures. The auto- 
mobile Is recognized as a major source of these primary pollutants In urban 
areas. However the control of automobile emissions has Increased the relative 
Importance of other sources of primary pollutants* such as airports. It Is 
now possible to estimate quantitatively the relative Importance of emissions 
from these other sources t^cause of our understanding of the photochemical pro- 
cesses that lead to oxidant formation. Perhaps more 1nq>ortant1y* this new 
knowledge makes It possible to assess the effects of the mixing of pollutants 
from a ntmiber of sources. Production or destruction of oxidant may be 
enhanced In an air parcel between two pollutant sources. Knowledge of 
such enhancement effects would be Important for consideration of control 
strategies proposed for either source or for consideration of the location 
of sources In an airshed area. 

The ratio of hydrocarbons to nitrogen oxides can be an Important factor 
in the production of ozone. Very high or very low ratios can lead to low 
ozone concentrations. At an Intermediate ratio much more ozone can be gener- 
ated from similar amounts of pollutants. In addition to differing from auto- 
mobile emissions In hydrocarbon c(xiX)os1t1on, jet aircraft emissions (especially 
from the Idle-taxi mode) have a higher hydrocarbon/NO^ ratio. This report 
shows that ozone enhancement effects are seen for the combined emissions of 
Jets and automobiles. 

The chemical behavior of jet emissions I'l the production of photo- 
chemical smog Is to be assessed In this study* Hoiimver* chemical effects 


t 
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•lone are exanlnedi the liopact on air quality of physical processes such as 
transport* dispersion* and dilution Mas not studied. The main elements of 
the effort vwre: 

> Estimation of emissions of Jet aircraft—:: literature survey 
of the types of compounds and their relative aixindance. 

> SlnNilatlon of the dynMics of photoch«n1cal smog— computer 
simulations using a ^nerallzed kinetic mechanism In static* 
contimwusly Irradiated air. 

> Evaluation and assessment of the sensitivity of the kinetic 
mechanism— assessment of the Impact on oxidant formation of 
variations In composition of miitted hydrocarbmss deter- 
mlnatim of sensitivity of prediction to variations In highly 
uncertain kinetic rate parmneters. 

> Estimation of ozone formation In air parcels containing mix- 
tures of auto and Jet Milsslons— conq^uter simulations that show 
enhanced ozone formation resulting from chemical Interaction 

In some cases. 


II EXHAUST EMISSIONS IN THE VICINITY 
OF AN AIRPORT 


Major fractions of the local hydrocarbons and nitrogen oxides observed In 
the vicinity of airports are attributable to emissions from aircraft (see Table 
1). Aircraft emit other compounds as well, such as water, carbon monoxide, and 
hydrogen, but at typical concentration levels these confounds do not affect the 
formation of ph..tochem1cal oxidant and are not considered In this study. 

Before the effects of aircraft emissions on local air quality can be 
examined the published data on aircraft exhaust emissions must be converted 
Into concentrations of pollutant species in air. The following data are 
required: 

> The rate of emission of total 'hydrocarbons, NO, and NOg from aircraft 
engines In the Idle, landing, and takeoff modes (the LTO cycle). 

> The time spent In each of the modes of operation by different 
aircraft at typical airports. 

> The fraction of the total hydrocarbons represented by each type of 
hydrocarbon (alkane, alkene, or aromatic). 

A. THE LTO CYCLE 

Aircraft operations In the Immediate vicinity of an airport can be 
separated into three modes: 

. * 

> Approach from the Inversion base to landing. 

> Landing, idling, and takeoff operations at the airport. 

> Climbout from the airport to the inversion base. 


Recent studies by Segal (1975) have shown that a typical worst case height 
of the Inversion base Is 150 meters. One kilometer Is the value assumed by 
the Environmental Protection Agency. 


Table 1 


PORTION OF POLLUTANTS ATTRIBUTABLE TO 
JET AIRCRAFT EMISSIONS 


Airport 

Hydrocarbons 

Nitrogen 

Oxides 

Atlanta. Georgia 

69 % 

78 % 

Los Angeles International 

72 

— 

O'Hare International. Chicago 

69.4 

74.7 


Reference 

Cirillo et al.(1975) 
LAAPCD(1971) 

Rote et a1.(1973) 


1 

O'Hare International Airport, one of the busiest airports In the world, 
serves as a central stopping point for passengers making transfers to other 
flights. Because transfers do not Involve automobiles, we expect that the 
percentage of pollutants due to automobiles should be low compared to other ■ 
airports. Rote et al. estimated the access vehicle percent at 10 percent. 


In order to prevent detonation from a cold running engine, aircraft use an 
enriched fuel-air mixture during Idle and taxi modes. The emissions from air- 
craft In these modes have high hydrocarbon concentrations and low nitrogen 
oxide (NO^) concentrations. During flight, combustion of a normal fuel-air 
mixture In a hot engine produces emissions having low hydrocarbon concentra- 
tions and high NO^ concentrations. Table 2 shows the emissions from different 
e.igines In various modes. Although one of the ways to qualitatively character- 
ize an operational nrade can be by a hydrocarbon/NO^^ ratio, the values of this 
ratio differ widely for different engines operating In the same mode. The 
ratio of total hydrocarbon emissions (THC) to NO^^ emissions In the taxi mode 
ranges from 8.0 for an Allison 501 -DB engine to 140.0 for a Pratt and Whitney 
JT-3D engine. A detailed discussion of these ratios will be presented In a 
later section. 

An examination of Table 2 reveals that for nx>st of the engines the 
largest amounts of hydrocarbon are emitted during the taxi -idle and approach 
modes. In Table 3, the typical time spent In each of the different opera- 
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Table 2 

CHISSION LEVELS FRON AIRCRAFT ENGINES 
IN VARIOUS OPERATIONAL MOOES 


Pollutant Species 


Enoine 

Mode 

HC 

(moles Carbon per enqine) 

— 

(moles NO.j oer engine! 

THC/N0„* 

Pratt A Whitney 

Taxi 

617.42 

4.16 

14S 

0T-3D 

Idle 

51.52 

0.39 

130 


Takeoff 

0.76 

8.12 

0.09 


Approach 

181.82 

19.3 

9.25 

Pratt & Whi tney 

Taxi 

67.05 

3.58 

18.4 

JT-80 

Idle 

5.68 

0.29 

19.2 


Takeoff 

2.65 

4.93 

0.53 


Approach 

12.12 

11.90 

1.00 

Pratt & Whitney 

Taxi 

46.97 

2.13 

21.6 

JT-12 

Idle 

3.79 

0.19 

19.6 


Takeoff 

0.38 

1.45 

0.26 


Approach 

1.14 

2.13 

0.53 

Allison 

Taxi 

36.36 

4.26 

8.4 

501- DB 

Idle 

2.65 

0.29 

8.97 


Takeoff 

3.03 

0.77 

3.87 


Approach 

11.74 

4.84 

2.39 

Pratt & Whitney 

Taxi 

245.45 

2.13 

113.12 

R-2800 

Idle 

16.29 

-- 

— 


Takeoff 

189.39 

— 

— 


Approach 

1091.29 

1.74 

615.72 

Continental 

Taxi 

18.94 

0.19 

97.86 

10-520-A 

Idle 

2.27 

" 

— 


Takeoff 

4.17 

-- 

— 


Approach 

25.76 

1.74 

14.5 

General Electric 

Taxi 

16.29 

0.68 

23.7 

CT-58 

Idle 

1.14 

0.10 

11.2 


Takeoff 

— 

— 

— 


Approach 

2.27 

4.74 

0.47 


♦Ratio of total caiton in moles to total moles of NO 2 . 
Source; Northern Research (1968) 






Table 3 

TYMCAL TIHE IN MODE FOR UNDIN6 TAKEOFF CYCLE 
AT A METROPOLITAN AIRPORT 


Time In mode— ml nutes 


Aircraft 

Taxi-Idle 

Takeoff 

Cllmbout 

Approach 

Taxi-Idle 

Jundio jet 

19.00 

0.70 

2.20 

4.00 

7.00 

Long range jet 

19.00 

0.70 

2.20 

4.00 

7.00 

Medium range jet 

19.00 

0.70 

2.20 

4.00 

7.00 

Air carrier 
turboprop 

19.00 

0.50 

2.50 

4.50 

7.00 

Business jet 

6.50 

0.40 

0.50 

1.60 

6.50 

General aviation 
turboprop 

19.00 

0.50 

2.50 

4.50 

7.00 

General aviation 
piston 

12.00 

0.30 

4.98 

6.00 

4.00 

Piston transport 

6.50 

0.60 

5.00 

4.60 

6.50 

Helicopter 

3.50 

0 

6.50 

6.50 

3.50 

Military transport 

19.00 

0.50 

2.50 

4.50 

7.00 

Military jet 

6.50 

0.40 

0.50 

1.60 

6.50 

Military piston 

6.50 

0.60 

5.00 

4.60 

6.50 


Source: U.S. Environmental Protection Agency (1973). 


tional modes is tabulated for different aircraft classes. The table 
shows that aircraft spend the most operating time at airports in the 
taxi-idle modes. The data presented in Tables 2 and 3 show that most 
hydrocarbons emitted from aircraft in the vicinity of an airport are from 
the taxi-idle modes. 


B. HYDROCARBON COMPOSITION OF AIRCRAFT EMISSIONS 

The rate of formation of oxidants in photochemical smog depends on the 
reactivity of the mix of hydrocarbons emitted by the various sources of 
pollutants. Whereas data on total hydrocarbon emissions from jet engines 
are plentiful, very little information is available on the quantities of 
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Individual hydrocarbons in the emissions. The most complete examination 
of hydrocarbon emissions of which we are aware is that of Conkle* Lackeyi 
and Miller (1975). They used a combustor test rig designed to simulate the 
operating conditions of the T-56 engine used on the C-130 transport. Dif- 
ferent operating conditions were simulated by varying the combustor Inlet 
temperature and pressure as shown in Table 4. At a constant inlet pressure, 
exhaust hydrocarbon compositions were measured for different types of 
fuels. Table 5 shows these measurements of hot exhaust gas. There can be 
substantial differences In exhaust emissions from use of the various fuels. 
The amount of paraffins emitted varied from 1.6 percent for JP8 fuel to 23 
percent for JP4 fuel, while aromatics varied from 20 percent fo JP4 fuel to 
52 percent for JP5 fuel. Taule 6 shows the variation of the percent of 
the total hydrocarbons In the exhaust under the various operating conditions 
for JP4 fuel. In most cases, che fraction of the species type varied 
by less than a factor of 2 over the range of operating conditions for 
which tests were made. Paraffins, olefins, aromatics, and aldehydes 
seem to be the major hydrocarbon constituents in the idle mode (correspond- 
ing to 15 psig). After adjustment to 100 percent, their respective per- 
centage values are 58, 18, 7, and 17 percent. 

Lozano et al . (1968) performed a study on the emissions from T-56-A7, 
J-57-19W, and TF-33-P5 engines. Their results are presented in Table 7. 

OP-4 fuel was used throughout their study. The hydrocarbon composition 
emitted in the idle mode from the T-56 engine contains approximately 62 
percent paraffins, 24 percent olefins, 9 percent aromatics, and 5 percent 
aldehydes. The largest difference between this composition and that re- 
ported by Conkle et al. is the percentage of aldehyde, which differs by a 
factor of 3. The hydrocarbon compos i ton reported for the exhaust emitted 
from the J-57 engine in the idle mode contains 49 percent paraffins, 24 
percent olefins, 25 percent aromatics, and 2 percent aldehydes. For the 
TF-33 engine, the corresponding values are 58 percent paraffins, 31 percent 
olefins, 8 percent aromatics, and 3 percent aldehydes. 
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Table 4 

SINGLE-COMBUSTOR RIG OPERATING CONDITIONS 


Nominal 
Power Setting 

Rig Pressure 
JrsM __ 

Inlet Temp. 
®C 

Fuel /Air 
Wt. Basis 

Fuel 

Preliminary 

75 

93 

0.0076 

JP4 

Moderate PR* 
Simulated idle 

33 

166 

0.0070 

JP4 

Low PR 

Simulated idle 

15 

93 

0.0073 

JP4 

Moderate PP 
Simulated idle 

33 

169 

0.0083 

JP5 

High PR 

Simulated idle 

50 

204 

0.0079 

JP5 

Moderate PR 
Simulated idle 

33 

164 

0.0071 

JP4 

Moderate PR 
Simulated idle 

33 

166 

0.0072 

JP8 

No fuel flow 
(Background) 

33 

— 

— 

— 


* PR--pressure ratio. 

Source: Conkle, Lackey, and Miller (1975). 
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Table 5 

DISTRIBUTION OF EXHAUST HYDROCARBONS BY LUMPED 
SPECIES FOR VARIOUS FUELS 


(at 33 PSIG Combustor Pressure) 


Species 

★ 

JP4 

JP4* 

JP5 

OPS 

No Fuel' 

Paraffins 

26% 

23% 

3.4 

1.6 

0 

Olefins 

24 

20 

25 

39 

0 

Diolefins 

3.2 

0.26 

t® 

0.45 

0 

Napthenes 

t 

2.3 

3.8 

0.45 

0 

Aromatics 

14 

7.2 

5.3 

2.3 

0 

Aldehydes 

20 

21 

52 

36 

0 

Alcohols 

1.9 

3.6 

0.57 

7.5 

82 

Ketones 

7 

7 

6.3 

10 

18 

Ethers 

1.6 

2.1 

3.4 

2 

0 

Esters 

0 

3.1 

t 

0 

0 

Nitrogen- 

containing 

compounds 

1.6 

0 

0.19 

t 

0 

Halogen- 

containing 

compounds 

1.6 

8.8 

0.19 

0.23 

0 

Lactones 

0 

0 

0.19 

0.23 


TOTAL HYDROCARBONS 

3. 17 

ppm 3.89 

ppm 5.27 ppm 

4.41 ppm 

0.11 ppm 


*Duplicate combustor conditions 
■•■Background sample (no fuel) 

^Trace, concentration less than 0.001 ppm 

Source: Conkle, Lackey, and Miller (1975) 
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Table 6 

DISTRIBUTION OF EXHAUST HYDROCARBONS BY LUMPED 
SPECIES AT VARIOUS INLET PRESSURES 


Combustor Inlet Pressure (pstq) 


Species 

15 

33 

33 

50 

75 

Peraffins 

SU 

23X 

26X 

31X 

27* 

Olefins 

16 

20 

24 

23 

47 

Diolefins 

0.2 

0.2 

3 

3 

t 

Napthenes 

2 

2 

t 

7 

2 

Aromatics 

6 

10 

14 

5 

14 

Aldehydes 

15 

21 

20 

23 

7 

Alcohols 

2 

4 

2 

0 

t 

Ketones 

3 

7 

7 

t 

t 

Ethers 

1 

2 

2 

3 

t 

Esters 

0 

3 

0 

0 

3 

Nitrogen- 

containing 

compounds 

0 

0 

2 

t 

0 

Halogen- 

containing 

compounds 

3 

9 

2 

4 

t 

Lactones 

0 

0 

0 

1 

0 

TOTAL HYDROCARBONS 

23.05 ppm 

3.87 ppm 

3.17 ppm 

0.74 ppm 

0.59 

Source: Conkle, Lackey> 

and Hiller (1975) 





Table 7 

POLLUTION EMISSIONS FROM JET AIRCRAFT 


‘ Power Setting and Fnqlne Tyne 

Take-off Cruise and approach Idle 


Pollutant 

T>5S 

J-57 

TF-33 

T-S6 

J-S7 

TF>33 

I:Si 

J-S7 

TFoa 

Oxygen (S) 

•• 

16.7 

17.1 

•» 

17.5 

18.0 

.. 

19.0 

19.6 

Carbon dioxide {%) 

4.1 

?.3 

2.7 

3.2 

1.5 

2.1 

2.4 

1.0 

0.9 

Carbon nionoxide (ppjn) 

34 

32 

7 

40 

55 

30 

109 

130 

195 

Oxides of nitrogen as 

K0> (ppn) 

NUnc ex tie (ppm) 

43 

59 

27 

27 

39 

15 

12 

13 

n 

37 

44 

25 

•• 

30* 

13 

— 

6 

9 

Total hydrocarbons 
(as C atO!T.s)(ppni) 

5.5, 

S 

7 

2.5 

5 

42 

101 

152 

700 

Olefins as C atoms 

(ppm) 

•• 

•• 

•• 

-- 

•• 

-• 

25 

38 

220 

Aromatics as C atoms 

(ppm) 

— 

— 

•• 

mm 

mm 

mm 

10 

39 

60 

Total aldehydes as 
HCKO (ppm) 

4.1 

0.8 

0.6 

2.0 

0.8 

0.3 

4.8 

2.5 

21 

Formaldehyde (ppm) 

1.1 

0.5 

mm 

1.9 

o.$ 

— 

3.5 

2.4 

mm 


Source: lo^ano et (1968) 
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In another study (Groth and Robertson, 1974), total and unreactive 
hydrocarbons were measured under various conditions for several Pratt and 
Whitney engines currently in operation on many commercial aircraft. The 
unreactive hydrocarbons consist of paraffins, which are relatively unroactive 
in smog-producing reactions with compared to olefins and aldehydes. 

The results of Groth and Robertson are shown in Tables 8(a) through 8(d). 

The unreactive hydrocarbons can reach as high as 30 percent during the 
idle mode (the idle mode is represented by the lower values of engine thrust). 
Groth and Robertson also determined the percentage of aliphatics, aromatics, 
and oxygenates in their study. Their results are summarized in Table 9. 

Since aliphatics contain olefins, alkynes, and diolefins, and oxygenates 
consist of aldehydes and ketones, we need tc know the percentage of alkynes, 
diolefins, and ketones emitted from these engines. For our study, we as- 
sumed the percentages are the same as those for a T-56 engine (Table 6). 

Using Tables 6 and 9, we calculated the composition of hydrocarbons emitted 
from a JT-9D engine during the idle mode as 51.3 percent paraffins, 15 per- 
cent olefins, 29.4 percent aromatics, and 1.93 percent aldehydes. For a 
JT-3D engine, the corresponding calculated values are 60.8 percent paraffins, 
17.8 percent olefins, 18.3 percent aromatics, and 0.8 percent aldehydes. 

The three studies presented in this section are the only studies we 
were able to find that contained extensive measurements of hydrocarbon 
compositions. Nevertheless a rather narrow range of hydrocarbon composition 
is indicated for emissions from different types of engines (e.g., turbojet, 
turbofan, and conventional jet), different operating conditions (e.g., tempera- 
ture, fuel-air ratio, and pressure), and different types of fuels (e.g., JP-4, 
JP-5, and JP-8). In Table 10, we sunmarize the hydrocarbon compositions from 
these three studies. 

C. AIRCRAFT EMISSIONS OF NITROGEN OXIDES 

Aircraft use a higher fuel-air ratio in the taxi-idle mode than in other 
modes, which leads to low combustion temperatures. Also the amount of nitro- 
gen oxides (NO^) emitted in the taxi-idle mode is low compared to the amount 






_ 
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Table 8 

HYDROCARBON EMISSIONS FROM JET ENGINES 
AT VARIOUS THRUSTS 


(t) JT'W tnsiM With JPS Fu*1 (b) JT*W tnglnt With Jtt A FmI 


Tlirwtt 

Totil Hydrocerbont 

Unreective Kydrccarboni 

Thrust 

Totil Hydrocifbont 

Unreeetive Hydrocarbons 

__ b-«vndi) 

(ppm) 

PPffl 

percent 

(pounds) 

(ppm) 

ppm 

percent 

199) (tubldle) 

82) 

)61 

4).7 

)9622 

2.0 

0 

0 

1626 (idle) 

407 

12S 

)0.8 

40)66 

1.6 

0 

0 

)419 

219 

68 

26.6 

42126 

1.6 

0 

0 

98)9 

89 

26 

28.1 

)7668 

2.0 

0 

0 

112)6 

60 

9 

16.0 

)49)4 

1.8 

0 

0 

19)47 

61 

7.4 

14.6 

16162 

).) 

0 

0 

2)16S 

44 

6.6 

14.8 

7100 

)6.6 

4.6 

12.) 

28626 

)6 

4.6 

12.8 

4000 

10). 6 

16.6 

16.0 





3040 

174 

)4.6 

19.8 


(c) 

JT-30 tnglnt With Riit Of Mtteriilt 







turner Cent And JPS Fuel 


(d) JT-3D Engine With Smokelett 

Sumer Cans And 0P6 Tue) 

Thrust 

Total Hydrocarbons 

Unreectivc Hydrocarbons 

Thrust 

Total Hydrocarbons 

Unreactive 

Hydrocarbons 

(pounds ) 

- 

PPffl 

percent 

(pounds) 

(ppm) 


percen’t 

760 (idle) 

666 

122 

18.4 

760 (Idle) 

477 

126 

26.3 

2100 

226 

63 

23.6 

2000 

110 

20 

18.2 

7S00 

)0 

8 

26.7 

• 7500 

8 

1 

12.6 

16800 

14 

0 

0 

16000 

6 

0 

0 

17900 

6 

0 

0 

18600 

6 

0 

0 

13000 

14 

0 

0 

13000 

2 

0 

0 • 

9960 

10 

0 

0 ' 

10000 

2 

0 

0 

3000 

8) 

8 

9.6 

3000 

32 

8 

26.0 

860 

665 

100 

16.3 

860 

410 

110 

26.8 


% 

i 

1 


Source: Groth end Robertson (1974) 
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TibU I 

liqulO CHROMT06RAPH ANALYSIS OP EXHAUST HYMOCADBONS 


CnOln. 

Twt 

Thrwtt 

tlbl 

ToUl 

Hydrocirboni (pipi) 

JT>M) 

2.702 (tub 

141.) 

ISS 

27.S0S 


S.4 

JT.JD 

2.700 (tub 

141.) 

170 

12.700 


30 

JT.30 

700 


1200 

3.SOO 


244 


A11ph«t1es 

■,.iBI£gnU 


71.1 

43.0 


10.7 

70.9 


OMyctnttti 

^ kiQ 1*^ {fitf 


aSl 


8SJCL. 

fittest! 

11 

39.4 

3.8 

2.S 

2.4 

44.4 

0.6 

11.1 

33 

19.4 

4.2 

2.S 

16 

63.3 

1,1 

3.7 

220 

18.3 

12.0 

1.0 

46 

18.9 

6.3 

2.2 


• Ptrctnt of total hydroearboni; assumption mada 
that 1 ppn ■ 1 111 1*» • 1 u9 !•* for coa^rlson 
purposas only. 

Sourta: Groth and Robtrtson (1974). 


Tibia 10 

COMPOSITION OF EXHAUST MYONOCAPBONS DUPING IDLE MODE 
(in pareent) 


En gine Typt 


Soeclis 


T-56-A7t 

0.57-TJWt 

TF-33-PSt 

UT>9DI 


Par# f fins 

68 

62 

49 

58 

51.3 

60.8 

01 af Ins 

18 

24 

24 

31 

15 

17.8 

Aromatics 

7 

9 

25 

8 

29.4 

18.3 

Aldehydes 

17 

6 

2 

3 

1.93 

0.8 

fuel used 

JP4 

JP4 

JP4 

JP4 

UP5 

UPS 


*Conk1i at i1. (1975) 
tLosano at i1. (1968) 

IGroth and Robartson (1974) 
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of NO^ emitted during the take-off and approach operational modes because of 
’ the power requirements Involved. Tables 2 and 7 show the amount of 
emitted In each of the different operational modes. The conA>1ned taxi and 
Idle mode: emit the lowest amount of NO^ for the following aircraft engines: 
Pratt and Whitney JT-30, Pratt and Whitney JT-8D, and General Electric CT-58. 
Similar amounts of NO^ are emitted In the different operational modes for the 
other aircraft classes. 

The results of the study by Lozano et a1. (1968) (Table 7) show that 
NO^ emissions are lowest In the Idle mode and highest during the take-off 
mode for the three engines studied. Table 11 presents a summary of emission 
factors for a number of different aircraft engines compiled by the Environ- 
mental Protection Agency (1973). Again* emissions are lowest during the 
taxi-idle mode. 

As with hydrocarbon compositions* very few studies have been made on 
the composition of the nitrogen oxides emitted from aircraft engines. The 
ratio of NO to NO^ varies between different types of engines. Lozano et al. 
(1968) reportc^d that roost of the emitted was In the form of nitric oxide 
(NO); they found that nitric oxide varied from 82 to 93 percent by volume 
of the total NO^^ emitted from a TP-33 engine (Table 7). For the J-57 en- 
gine, NO varied from 62 to 76 per{:<»nt of total NO^ depending on the opera- 
tional mode. Lo/ano et al. also found the percent composition of NO was 
greatest In the take-off mode and Irwest In the Idle mode. 

In a detailed study by Bogdan and McAdams (1971)* aircraft exhaust 
emission measurements performed by different organizations are summarized. NO 
and NO 2 coticentrations were measured Individually for many of the test engines. 
The value of th: NO/NO 2 ratio was found to ran^ from approximately 1 to 20 for 
most engines. Table 12 snevs the NO^ emissions for some engines; the average 
NO/NO 2 ratio IS approximately 9.7. 


Table 11 

MODAL EMISSION FACTORS 
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btLntfjL 

C.rlMii 


(NO, 4, NO,) 

S,IM 

i«f.Ui«tiSSL 


.jm 

mt 

I^Of 0«i/fir 

lo/>y 

Jozoc 

-Mic 

jMOc. JaOt 

Prt%% 1 UMtiit/ 

jH) 

Taxl-ftfit 

}.m 

701 

m.o 

40,3 

37.3 

13,4 

0,00 

3.7$ 

3,3 

1.0 


I7.0SI 

7.73$ 

1,30 

3.70 

3,0$ 

1.34 

730,0 

X7.0 

3,7$ 

1,7 

ft 

14417 

$.404 

11,7 

$.31 

3.0$ 

1,30 

410.0 

300.0 

4,0 

1,0 

AK'^-iCfc 

S.I04 

f.3$l 

X.O 

14,0 

3,00 

1.M 

$4,1 

34.1 

3.3 

1,0 

ttutfAl Clac- 
trfc m 

j«t) 

TixUUli 

l»W 

407 

01,7 

33,0 

10.4 

7.0 

3.0 

1.03 

0.04 

0.03 

T«ki«rr 

I3»44t 

0.1CO 

0,7 

3,04 

1,3 

o.$o 

$40,0 

34$,0 

0,04 

0.34 

Cl 1169^^ 

1M00 

$.171 

0,0 

3,00 

1.3 

0,$0 

333,0 

1$1,0 

0.14 

0.34 


<»K4 

1.014 

10,0 

0.44 

1,0 

0.00 

173.0 

70.$ 

0.44 

0.30 

^ilt 1 NMtiiiQr 
JT*I9 

(iMif rif^ Jtl) 
TaiMtfIt 

in 

300 

100,0 

40.4 

00,0 

44.7 

1.43 

0.040 

0.4$ 

0.30 

Tikaoff 

10,13$ 

4.01$ 

13,3 

$.00 

4.0$ 

3.11 

140.0 

07.1 

0,» 

3.7 

€ll«tettt 

••$$$ 

4.003 

IS, 3 

0,04 

4,03 

3.33 

M.3 

43,0 

0.1 

3.0 


4«13l 

1,077 

30.7 

10.0 

7.04 

3,H 

31.0 

0.00 

0,0 

3.0 

Pr%ii 1 
4T-X 

(L09f r«AH Jtt) 
TaxMtflt 


$43 

03.0 

43.0 

03,3 

41.0 

3.40 

1.13 

0.40 

0.10 

Tiltoff 

10.113 

4,010 

0,04 

4.10 

0.0$$ 

om 

110.0 

$4.0 

O.M 

3.0 

CIMoitC 

••sot 

3,000 

10,0 

7,30 

0.001 

0.4M 

04.7 

X.4 

0.3$ 

3,0 

ApymcH 

4.11$ 

1,007 

40.0 

33,3 

0,30 

3.7$ 

33,3 

10.$ 

3,» 

1.1 

1 tAiltMir 

^*4A 

(t0A9 rAft9* j«t) 
T«xM4l« 

1.3tf 

030 

03,0 

’30,$ 

04.0 

30,4 

3.71 

1.33 

1.3 

0,$4 

TtMoff 

is.sn 

7.030 

10.0 

0.S3 

0.074 

0.300 

330,0 

107.0 

21,0 

0.1 

Clliibeiit 

13.0$$ 

$.037 

10.3 

0.30 

1.37 

0.$70 

1$S,0 

70.3 

30.0 

0.1 

Aa^fMcA 

$.H4 

3,710 

20,3 

11.0 

3.03 

1.74 

3$.0 

10.3 

0.0 

3.7 

(efieril tUc- 
trU C49M 

UoMf riRfe jft) 
TiiM4U 

1.001 

4$4 

03.0 

:: 0 

37,3 

13,4 

1,17 

0,713 

1.3 

0,$0 

UUafr 

t.OlO 

4.110 

30.1 

13 ; 

O.SSO 

0.3S3 

111.0 

$0.3 

1S.0 

0,0 


o.m 

3.760 

30,0 

•j.l 

0.SO3 

0.304 

74.0 

33.0 

u.o 

0.0 

ApafMcii 

3.717 

1.713 

43,0 

10.4 

3.43 

1,10 

17,0 

0,07 

$.0 

3.3 

fnxi i imUiiejf 
j«l) 

laaMilt 

f$l 

43$ 

33.4 

1$,3 

0.00 

3.71 

3,01 

1.33 

0.30 

0.10 

Ukfff 

0.7$$ 

3.071 

7,40 

1.40 

0.770 

0.3IJ 

100.0 

00.0 

3,7 

1,7 

tU9tm% 

7.317 

3.330 

O.H 

4.01 

0.031 

0.410 

131,0 

$0.0 

3,0 

1.3 


3.401 

1,140 

10.3 

0.30 

1,7$ 

0.704 

30.0 

14.0 

1.1 

0.00 
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Table 11 (Concluded) 



. fi«l rttf 

CmOw 
■W unl,. 


M«n 
(W, M M,) 

totM 

•.ntMUttf 

Am^m md ^ 


JtdS. 

Jl£S£. 


.jm 

IfcTOf^ l8/lf 


S^wtlU 
^.riHt Jtl) 
TotHdlf 

m 

100 

80.1 

17.3 

88.1 

30.0 

0.848 

o.m 

0.17 

8.877 

Tilttff 

Mil 

1*411 

14.1 

1.44 

8»8 

048 

1S3.0 

88.4 

18.0 

7.3 

aiiiiM 

i»iH 

1,811 

U.l 

8.84 

0.141 

0.110 

11S.0 

41.7 

10.0 

4.1 

^TMCl 

t«0SI 

1,184 

18.1 

17.7 

4.77 

1.81 

30.4 

13.0 

1.8 

8.88 

AiiitMi m-Aii 

(Air cirritr Uflo* 
pfOfi nil. trm- 
Hrtl 
TuMdlt 

Ail 

tH 

8.74 

3.81 

7.38 

3.IS 

1.73 

0.180 

1.8 

8.73 

Tiktirf 

t.m 

1,081 

1.77 

1.71 

0.440 

0.700 

17.8 

11.7 

3.7 

1.7 

OlMloyt 

t.m 

m 

1.40 

1.S4 

0.3H 

0.181 

17.7 

10.1 

3.0 

1.4 

ApiTMili 

M4I 

IfO 

1.48 

1.18 

0.378 

0.148 

7.37 

3.37 

3.0 

1.4 

AIHim TM-A7 

^If cirrltr tuflo* 
tfopi nil. triii* 
P9tt) 

iMl-IAIt 

m 

141 

11.1 

8.84 

8.47 

t.n 

M8 

0.880 

1.8 

0.73 

tlkttff 


841 

MS 

0.871 

0.4^ 

0.18S 

17.8 

10.4 

3.7 

1.7 

ClInlM 

MOS 

US 

1.01 

1.37 

0.478 

0.718 

71.7 

8.87 

3.0 

1.4 

AfpfMCl 

t.iu 
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1.17 

1.18 

0.S17 

0.7H 

7.78 

3.83 

3.0 

1.4 

A<r«i»iftl 

W-Ml** 

(fin. iVlitiM 
tdftofrvp) 

TixMdli 

Ul 

HA 

I.SI 

• 
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0.878 

0.3H 

0.8SS 
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m 
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0.178 

O.MS 

0.07S 
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1.8S 
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O.OS3 
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m 
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2.S8 
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0.787 

0,8 
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Ttlf4jw/C©«- 
tififbti) o-ioo 
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>.$7 
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m 
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NA 

8A 
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0.377 
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KA 
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pliUfi) 
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11.0 

l.fO 

11.1 
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NA 

Tike^ff 
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70,8 

37.7 

1.48 
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M 

NA 
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78.8 

1.31 

0.584 

0.37S 
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m 

NA 
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13.1 

10. S 

^4<3 

11.0 
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0.011 
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AfUlyiU ef Aircraft C^beuit fowll lefioret^f/, )rtc. W,Y. 

frwred fnr t»»* T;^4MrUi >#rk. U,C*t uftCtf C^trect 6l»ie<»ef 

^04-^n. (K-UttrW\, 
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All i<f itffltrs (chedyled for cy«%fr4len of !• bfaiekeltift If 
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evetUbIt 


Enpln* 

Pr«tt and Ubitney JT-30 

rratt and Whitney JT-80-1 

Pratt and Whitney JT-9D 

Pratt and Whitney J57-P10 

Pratt and Whitney J-79 
Pratt and Whitney T-56-A7B 


Table 12 

EMISSIONS FOR SOME ENGINES 


Thrust (lbs) 

NO (ppnV) 

NO? (ppmV) 

Ratio NO/NO? 

1777 

7.2 * 

8.6 

0.8 

€778 

27.9 

6.9 

4.04 

14628 

80.5 

2.9 

27.8 

9660 

61.0 

5.0 

12.2 

11480 

84.0 

5.0 

16.8 

800 

4.0 

7.0 

0.57 

2482 

9.6 ■ 

11.1 

0.9 

3493 

13.8 

15.3 

0.9 

10968 

62.4 

13.0 

4.8 

430 

7 

2.9 

2.4 

6080 

•30 

3 

10 

7940 

43 

3 

14.3 

9420 

67 

3 

22.3 

275 

8.0 

3.0 

2.8 

6650 

39.0 

4.0 

9.8 

9800 

62.0 

6.0 

10.3 

2237 

71.0 

7.0 

10.1 

3114 

103.0 

5.0 

20.6 

3443 

109.0 

5.0 

21.8 


Source: Bogdan and McAdams (1971) 
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D. HYDROCARBON/NOj^ RATIOS FROM AIRCRAFT EMISSIONS AND AUTOMOBILE EMISSIONS 

Con?)uter simulations of the effects of mixing aircraft and automobile 
emissions require hydrocarbon/NO^ ratios for emissions from aircraft and auto- 
mobiles. The results of these simulations are presented in Chapter V, but the 
emissions and HC/NO^ ratios are discussed here for easy comparison with the 
previous section. 

In order to tabulate a consistent set of hydrocarbon/NO^ ratios from 
available data, the following conditions were imposed: 

> All data were converted to moles of carbon for hydrocarbons and 
moles of NO (or NO 2 ) for NO^. 

> Ratios of hydrocarbon to NO^^ were required for the whole LTO 
cycle. 

Conversion factors are discussed in the Appendix. 

1. Hydrocarbon/NO^ Ratios from Aircraft Emissions 

Due to the broad range of values for hydrocarbons and NO^ emitted from 
different aircraft engines using different types of fuel and different opera- 
tional modes, a comprehensive analysis of reported data was performed. Table 
2 presents hydrocarbon/NOj^ ratios tabulated for each mode of the LTO cycle 
from data obtained from Northern Research (1968). Table 13 lists the total 
hydrocarbon and NO^ emitted during the LTO cycle. These values were obtained 
by combining the emissions from the different modes. All emissions were 
multiplied by a time factor taken from Table 3. The time factors were obtained 
by calculating the relative percentage of time spent in each mode. For the 
seven classes of engines, the hydrocarbon/NO^ ratios range from 'v- 400 to 3.6. 
For the three engines studied by Lozano et al. (1968), hydrocarbon/NO ratios 

A 

varied from 2.7 for the T-56 engine to 27.0 for the TF-33 engine (Table 14). 
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Table 13 . 

SUMMARY OF ENGINE EMISSION LEVELS DURING AN LTO CYCLE 


Engine 

HC 

(moles of carbon 
per engine) 

NOx 

(moles of NO 2 
per engine) 

THC/NO; 

Pratt and Whitney 
0T-3D 

851.5 

32.0 

26.1 

Pratt and Whitney 
JT-8D 

87.5 

20.7 

4.1 

Pratt and Whitney 
JT-12 

52.3 

5.9 

8.7 

Allison 
501 -DB 

53.79 

10.2 

5.2 

Pratt and Whitney 
R-2800 

1542.4 

3.87 

391.0 

Continental 
10- 520- A 

51.1 

0.39 

129.0 

General Electric 
CT-58 

19.7 

5.5 

3.5 


* 

Units of moles of carbon/mol es of NO 2 
Source: Northern Research (1968) 

Table 14 

POLLUTION EMISSIONS FROM JET AIRCRAFT 
DURING AN LTO CYCLE 


Engine 

HC (ppm as C) 

NOx (ppm as NO 2 ) 

THC/NO 

T-56 

58.2 

21.2 

2.7 

J-57 

87.3 

27.6 

3.14 

TF-33 

404.9 

14.7 

27.0 


Source: Lozano et al. (1968) 
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EPA (1973) has reported modal emission factors for 13 different aircraft 
engines with uses ranging from jumbo jets to air carriers (Table 11). Table 
15 shows the emissions during the entire LTO cycle and the hydrocarbon/NO 

A 

ratio for each engine. The values for the hydrocarbon/NO^ ratio range from 3.5 
to 70.0 (moles as C/moles as NO 2 ). 

Broderick et al. (1971) performed a survey of aircraft emissions monitor- 
ing requireirents. Table 16 lists the types of jets and engines considered In 
their study. Table 17 shows the emission Index of the jets listed In Table 16. 
These values were cited by Broderick et al. from a study by Northern Research 
(1968). Again, there Is a broad range of values for the hydrocarbon/NO^ 
ratios— from 33.0 to 137 (units of iw)les of carbon per moles of NOg). 

Naugle (1974) surveyed the measurements of exhaust emissions from 
military aircraft. Table 18 summarizes some of the aircraft emissions data 
tabulated by Naugle. This table lists the emissions during the entire LTO 
cycle. The hydrocarbon/NO^ ratios range from 2.4 to 22.4 for the four engines 
studied. Table 19 lists emissions per touch-and-go cycle. The touch-and-go 
operations are used as training methods at many Air Force bases. According to 
Naugle, these operations can give emissions "nearly as great as from landings 
and takeoffs at some airbases." The hydrocarbon/NO ratios are low compared 

A 

with the ratios in Table 18. 

Rote et al. (1973) reported emissions of aircraft at O'Hare International 
Airport. Table 20 lists their observations. The hydrocarbon/NO^^ ratios range 
from 4.0 for the superjets to 76.0 for long range jets. 

The six studies cited in this section report a broad range of values for 
the hydrocarbon/NO ratios. Therefore, for the computer simulations, an 

A 

averaged hydrocarbon/NO ratio was used. The averaged value for the hydro- 

A 

carbon/NO ratio is approximately 41.4. 

A 

The hydrocarbon/NO ratios tabulated for emissions during the LTO cycle 

A 

provide one measure of the hydrocarbon/ NO^^ ratio, but the gross amounts of 
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Table 15 

EPA MODAL EMISSIONS FACTORS DURING 
LTO CYCLE 

Ratio of Tota’ Hydro 


Engine 

Hydrocarbons 
(kg hr“l) • 

Pratt and Whitney JT-90 

10.7 

General Electric DF-6 

6.07 

Pratt and Whitney JT-3D 

38.5 

Pratt and Whitney JT-3C 

36.0 

Pratt and Whitney JT-4A 

25.2 

General Electric CJ-805 

10.7 

Pratt and Whitney JT-8D 

3.26 

Rolls Royce Sprey MK 511 

25.8 

Allison T56-A15 

2.84 

Allison T56-A7 

2.5 

Ai research TPE-331 

0.35 

Teledyne/Continental 0-200 

0.12 

Lycoming 0.320 

0.18 


NOv (as NO2) carbons to HOx 

(kg hr*^) (ppm as C/ppr as NO 2 ) 


12.1 

3.4 

16.5 

1.41 

3.18 

46.5 

3.41 

40.6 

5.30 

18.3 

2.66 

IS. 5 

4.74 

2.64 

3.51 

28.3 

1.22 

8.95 

1.56 

6.17 

0.51 

2.64 

0.01 

46.1 

0.01 

69.2 


Source: EPA (1973). 


Table 16 


AIRCRAFT CONSIDERED IN BRODERICK'S STUDY 


Type 

Long-Range Jet 
Medium-Range Jet 
Turboprop 
Piston 
Super jet 


Example 

707, DCS 
727. 737, DC9 
El ectra 

DCS, Corvair 440 
747 


Engine 

Pratt & Whitney JT-3D 
Pratt & Whitney JT-8D 
Allison 501 -DB 
Rolls Royce R2800 
Pratt & Whitney JT-9D 


Table 17 

EMISSION INDEX CITED FROM BRODERICK'S STUDY 


. _ Iype _ . 

Fuel Usage 
(lb hr"^l 

HC 

(1000 lb fuel) 

Long-Range Jet 

1090 

54.6 

Medium-Range Jet 

920 

12.7 

Turboprop 

587 

9.7 

Super Jet 


13.4 


Ratio of Hydrocarbons 
NOx to NOx 

(1000 lb fuel) (ppm as C/ppm as NO^) 


1.56 

134 

1.46 

33.5 

2.18 

17.1 

1.55 

33.3 


Table 18 


EMISSIONS PER LTO CYCLE 


Aircraft 

Engine 

Hydrocarbon 
(lbs /cycle) 

NO^ 

(lbs /cycle) 

HC/NO 

(moles as C/mofes as NO^) 

F-4 

J-79 

9.8 

14 

2.7 

F-104 

J-79 

4.4 

7.3 

2.4 

T-37 

J-85 

19 

3.7 

19.7 

T-38 

J-85 

25 

4.3 

22.4 

B-52H 

TF-33 

680 

170 

15.4 

KC-135A 

J-57 

198 

61 

12.5 


Source: 

Naugle (1974) 


Table 19 



• 

EMISSIONS PER TOUCH-AND-GO CYCLE 


Ai rcraft 

Engine 

Hydrocarbon 

(Ibs/cycle) 

NO 

(lbs /circle) (moles 

HC/NO 

as C/moTes as NOo) 

F-4 

J-79 

0.22 

5 

0.17 

F-104 

J-79 

0.15 

3 

0.20 

T-37 

J-85 

00 

• 

in 

2.8 

7,96 

T-38 

J-85 

4.6 

2.4 

7.37 

B-52H 

TF-33 

6 

52 

0.44 

KC-135A 

J-57 

6 

25 

0.92 


Source: Naugle (1974) 
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Table 20 

EMISSION RATES PER LTO CYCLE 


Aircraft 

Engine 

Hydrocarbons 
lb hr-1 engine"! 

NOx 

lb hr"l engine"! 

HC/NOj^ 

(moles as C/ 
moles as NO 2 ) 

Jumbo 

JT-9D 

32.7 

30.7 

4.1 

Long Range 

JT-3D 

144.6 

7.3 

76.2 

Medium Range 

JT-8D 

15.4 

8.7 

6.8 


CJ-805-3A 

46.2 

7.9 

22. S 

Short Range 

A-501-013 

10.3 

4.9 

8.1 


SPEY 511 

107.7 

8.7 

47.6 


Source: Rote et al. (1973) 


Hydrocarbons and NO^ emitted by all aircraft at an airport can provide another 
measure. Northern Research (1968) reported the total emissions from all air- 
craft operating at FAA-contro11ed fields as follows: hydrocarbons— 415 tons/ 
day; NO —42.8 tons/day. The hydrocarbon/NO^ ratio for these values Is 37.3 
(moles as C/moles as N0»). Rote et a1. (1973) reported emissions of 5.074 x 
lO"* tons/year for hydrocarbons (moles as C/moles as NOg) and 0.2931 x 10 
tons/year for NO , for a hydrocarbon/NO^ ratio of 66.6. LAAPCD (1971) 
reported the following emissions at the Los Angeles International Airport: 
10,685 tons/year of hydrocarbons and 1105 tons/year of NO^^. The hydrocarbon/ 
NOj^ ratio Is 37.2 (moles as C/moles as NOg). Averaging the ratio found for the 
Individual engines with the three ratios of gross emissions yields a value of 
45.7 (In units of moles as carbon/moles of NO^)* 

2. Hydrocarbon/NO Ratios For Automobile Exhaust Emissions 

Unlike aircraft engine operations during the taxi-idle mode, automobile 
engines (from 1972 to the present) run on a "lean" fuel-air mixture. Automo- 
bile exhaust emissions have low hydrocarbon concentrations and high NO con- 

A 

centrations. This is the opposite of aircraft emissions. Table 21 presents 
the results of a Bureau of Mines study (1973) on ten 1970-1973 automobiles 
using a 40 percent aromatic fuel. The HC/NO ratios are low compared with the 

A 

ratios for jet aircraft (Table 10). The EPA studies (Table 22) show slightly 
higher hydrocarbon/NO ratios for average emissions in a given year. An 

A 

Important feature shown In Table 22 Is that greater amounts of hydrocarbons 
are emitted at the slower speed than at the higher speed. 

Of the ten automobiles studied by the Bureau of Mines (Table 21), only 
the 1971 Vega shows a very high hydrocarbon/NOj^ ratio. Excluding the Vega, 
the range of values Is very narrow compared with the wide range of values 
found for aircraft emissions. The average hydrocarbon/NO ratio for automobile 

A 

emissions from these cars is approximately 2.2 (excluding the Vega). The 
higher ratios reported by the EPA (Table 22) have an average hydrocarbon/NO 

A 

ratio of approximately 4. 
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Table 21 

EMISSIONS FROM 1970 - 1973 AUTOMOBILES 


Vehicle 

Total Hydrocarbons 
(grams/fflile) 

N0„ 

(grams/fflile) 

HC/NOJ. 

1972 Oldsmoblle (car 151) 

1.79 

5.01 

^ A 

1.4 

1971 Ford Galaxle (car 707) 

4.88 

7.34 

2.6 

1971 Plymouth Fury (car 76) 

3.22 

8.52 

1.45 

1972 Ford Torino (car 769) 

2.16 

4.32 

1.92 

1970 Chevrolet Impala (car 595) 

3,58 

4.37 

3.15 

1970 Pontiac (car 400) 

5.86 

9.36 

2.44 

1970 Volkswagen (car 365} 

2.44 

6.09 

1.54 

1971 Vega (car 68) 

4.92 

1.88 

10.1 

1973 Ford Torino (car 146) 

2.79 

4.17 

2.57 

1973 Impala (car 110) 

2.16 

2.66 

3.12 


Note: Fuel used was 40 percent aromatic. 
^Moles C per mile / moles NO^ per mile. 

Source: Bureau of Nines (1973) 


Table 22 

AUTOMOBILE EMISSIONS BY YEAR AMD NATIONAL AVERAGE EMISSIONS 


Vehicle 
1971 Models* 

1971 (California)* 

1972 Models* 
National Average * 

a 45 mph 
9 19.6 mph 


Total Hydrocarbons 
(qrams/mllel 

7.2 

2.9 

6.6 

4.7 

14.6 


NO (as NO-) 
(grams /mile) 

5.4 


3.5 

5.4 


8.0 

4.6 


• Source: LW u9?3) 

Source: EPA (1975) 

^Holes C per mile / moles NO^ per mile. 


HC/NOj^ * 

5.14 

3.19 

4.7 


2.3 

12.2 


Ill COMPUTER SIMULATIONS OF PHOTOCHEMICAL SMOG 


An important part of this study is the prediction of oxidant forma- 
tion given an initial set of conditions. Ozone is considered to be the 
principal component of photochemical oxidants. The nwthoci for the mea- 
surement of oxidants is specific for ozone (EPA, 1971). Therefore, this 
study was confined to the ozone behavior predicted by the kinetic mechanism. 
Predictions were made by numerical integration of a set of differential 
rate equations that describe the kinetics of photochemical smog. The 
numerical integrations were performed using the CHEMK program developed 
by 6. Z. Whitten of Systems Applications, Incorporated on the CDC 7600 
computer located at the Lawrence Berkeley Laboratories. 

A. THE GENERALIZED KINETIC MECHANISM 

The kinetic mechanism used as a starting point for the computer simu- 
lations In this study Is an updated version of the generalized 39-step 
mechanism developed by Hecht, Seinfeld, and Dodge (1974). The mechanism 
Is presented In Table 23. This kinetic mechanism has been validated by 
Hecht et al. (1974a, b) using data collected in smog charrber experiments 
performed at the University of California at Riverside (UCR). The ©mis- 
sions data surveyed in the present study were used as starting conditions, 
as In a smog chanter experiment, for the computer simulations. The com- 
puter program (CHEMK) Is capable of simulating up to 200 rate expressions 
containing up to 50 species. The CHEMK program has been previously 
employed in the validation of kinetic mechanisms by Durbin et al. (1975) 
using UCR Smog chamber data. Simulating using the emissions data as If In 
a smog chamber experiment can be considered as an illustrative situation, 
because the system remains static, I.e., there Is no loss or gain of pro- 
ducts or reactants due to dilution nr replenishment of the system, and 


Table 23 

THE SH06 MECHANISM 
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— ■ 

NOj^t Oj Inorytnlc clMUtry 

0*0^«N«0, «N 

0, ♦'no* NOj ♦ 0| 

0 « NO « 

0, ♦ NOj NOj ♦ Oj 

MOj « NO - 

NOj « NO^ •> NjOj 
NjOj « N ■* NOj « KOj 
OK « • N»0, 

ON « NO « MIO^ 

NOj * NO * UOj ♦ OH 

HO^'NOj reconbinitton 

KO^ ♦ HOj •• HjOj 


N*tt CmtMt 
' <«*»'* »*«»•*) 

t.OI M 10** OPN** rtK** 
ts.t 

1.34 s 10* 

O.OS 

1.3 K 10* 

1.0 X 10* 
f.4 X 10** 

1.0 X 10* 

3.0 X 10* 

0.0 X 10* 


0.0 X 10* 


Or$8n1c oxfdAtlon reictlpns 

Alkxn* ♦ OK - Wj « MjO 3.4 x 10* 

AUXM 4 0 AO; « OH ** . 

AUene 4 OH ■* Ald»hyd» 4 ,hO; ♦ (1 - a)W; t.S * 10* 

AUen# 4 0, ■» OH 4 Aldehyd* 4 »RCOj 4 (| . »)HOj 4 (1 . «)C0 0.03 

AUent 4 0 - AO; 4 yWj ♦ tCO ♦ (1 - t)RCOj 5.3 x 10* 

Aroaxtlc 4 OH . M; 4 H;0 0.0 X 10* 

Aro«itk ^ 0 -• RO^ ♦ OH 120 

AldcKydt 4 OH - «H0; 4 (1 - 3.1 x 10* 


XftctloM of luiapod xoeclot 

RO; 4 NO * Aldehyde 4 0.5(1 ♦ t)HO; ♦ 0,5(1 ♦ i)CO 4 0.5(1 - c)RO; 4 
NCO, 4 NO « RO; 4 CO; 4 NO; 

NCO; 4 l«; .. CAN 
RAN » RO; 4 CO; 4 NO, 

RO; 4 HO; 4 RO;H 4 0; 

WO, 4 MO; 4 RtO,H 4 0; 


1.0 X 10* 

1.0 X 10* 

3.0 X 10* 

3.0 X 10** xilo"’ 
3.0 X 10* 

3.0 X 10* 


Table 23 (Concluded) 




lUtt CoAtUnt 


to Midttlon 

CO * ON ■* COj 4 HOj tSO 


fAotocNmlttry 
N 0 | « Nv * NO * 0 

KgOf * ho - ton 

Aldthydt 4 hv " fro<>wcti 

AiM/dt 4 • (1 * 4 0 . 4}K)^ 4 CO 

MNOf 4 Nv * ON 4 NO 

St/rftct rttctiofit* 

NO 4 NOj 4 HjO 4 ?HN0j 

• 1(0 4 NOj 4 HjO 
NjOj 4 H^O * ?hNOj 

Oj • Products 


k| Kin*' 

O.ON k, din*’ 

0 . 004 S k| Mln*^ 

0 . 004 k k| nfn*^ 

o.o; k, 

C.C « »0*^ pp»‘* H«*’ 
0.?0 
S N 10’* 

0.001 n(n’’ 


• • splU PO *• * AlMydt «nd RO ^ 0 ^ ** HO^ ^ • O.R» 

$ • ip) It betk^n pithiMiyi for e*p Mrjdfci) dtcomposU^on • 0 . 33 * 

Y « ipHt t^Ptwpen Inttrni) iftd Urr,in«) 0 idd^tfon • 0 ,S. 

I • ffictlrn cf itdehydts thit art fcrmaldt^df * 0 *S. 
c • fraction of RO^ that U 

*Thest depend on the turface/vo 1 u<nt and the niturt of the lurfactp 

liote: i and e are erpirka) “clotwre parameten.* The prelent values were 
determined from i 1 :?^^Ut 1 cnf of UCR data and the data themselves. 

f, and Y are functions of olefin. Their values «kere asslsned on 
the hails cf propylene values. 


Smxtf: 

Kiipfpfi and darvlA IW%) 
AtMnton Oi^ RUU (mS) 
filvit (197S) 

DditiA ft 0l. (137$) 

HfCht ft f1. (1374b) 
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solar radiation 1$ held constant for a 12-hour period. The hydrocarbons 
used In the validations of smog chanter data are butane* propylene* and 
toluene* which represent their respective general groups (e.g.* butane 
represents the alkanes In the reactions involving alkanes) for our Initial 
simulations. The closure parameters* a* a* y* d* and c are determined by 
the specific hydrocarbons used In the mechanism. For exanq>1e* the closure 
parameter a represents the fraction of the aldehyde produced which Is 
formaldehyde. If the alkane and alkene molecules are longer chain mole- 
cules, then the fraction of formaldehyde produced Is expected to be smaller. 
The value of S was determined from smog chamber data on a propy1ene/N0^ 
system obtained at the University of California at Riverside. The param- 
eters a* 8* and Y are functions of the alkene. 


The parameter a Is determined by the split between the unimolecular 
decomposition of alkoxyl radicals (RO*) and the oxidation of alkoxyl radi- 
cals. Although the Individual rates of these reactions are of little 
Importance In the kinetic mechanism* their ratio, a, Is Important In deter 
mining the course of the overall reaction. The parameter 8 Is determined 
by the pathways for the oxyperoxy b1 radical decomposition. The mechanism 
for this decomposition* developed by O'Neal and Blumstein (1973), Is now 
given (assuming the alkene Is propylene): 

Oj ♦ CH^-CHCHjV, CHjCH— CH2 (a moloronlde) . (1) 




) 1 


CH3CH 


CH, 



(a biradical) 


( 2 ) 


(3a) 


Blradlcal 


A 



(3b) 
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Whether reaction (3a) or (3b) occurs depends on the nature of the bi radical 
Intermediate. The a-keto hydroperoxide formed in Reaction (3) is in an 
excited state. Its fractionation leads to the following overall reaction: 


(4) 

(5) 

The split between Reactions (4) and (5) is determined by which hydroperoxide 
is formed in Reaction (2). Durbin et al. (1975) used a value of 0.33 for B 
to fit smog chamber data from University of California at Riverside. The 
parameter y is the fraction of carbons attached to the double bond in a 
monoolefin that are not terminal carbons on the chain. Consider the reactions 
of oxygen atoms with propylene and 2 -butene: 


Og + CH2=CHCH3 



0 

II 

H^CO + CH 3 CO 2 + OH* 


CH 3 CHO + CO + HO^ + OH* 


On 


* 



CH3CH^ CHO* 


0 

II 

CH3C* + CH^ 


> 


0 + CH3CH-CHCH3 


r 0 


CH3CCH2CH3 


CH3C- + CH3CH^ 


If we assume that the oxygen atom will react with equal probability with 
either of the carbon atoms in the double bond, y = 0.5 for propylene and 
Y = 1 .0 for 2-butene. 
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For the simulations performed In this chapter, we assumed a value of 
the NO^ photolysis rate constant of 0.25 ppm“^ mln”^ (an average value 
of for a typical solar spectrum). Some of the surface reactions listed 
In Table 23 have been shown to be very sensitive to changes in reactant 
concentrations and changes In rate constants (Durbin et al., 1975). These 
surface (or heterogeneous) reactions will be discussed In a later section. 

B. RESULTS OF THE COMPUTER SIMULATIONS 

A necessary input to the kinetic mechanism is a set of initial con- 
centrations of various groups of hydrocarbons. Four mixtures with differ- 
ent relative concentrations of alkanes, alkenes, and aromatics were 
constructed to represent the range of hydrocarbon compositions in jet 
exhaust emissions (Table 24). One of these mixtures and a value for the 
total hydrocarbon concentration were used in each computer simulation. 

The initial concentrations of hydrocarbons and used in the simulations 
are shown in Figure 1. For each mixture in Table 24, computer simulations 
were performed using each point in Figure 1 as an initial condition. 

Ozone profiles as a function of time were computed for each point. Iso- 
pleths of ozone concentrations were then constructed as a function of the 
hydrocarbon and concentrations at intervals of 1 , 2, 4, 8, and 12 hours. 
Due to the limited number of simulations, some interpolation was necessary 
in constructing the isopleths. 

Since the computer program is coded to accept initial concentrations (in 
ppm) of specific molecules, hydrocarbon concentrations reported as ppm of 
carbon had to be converted before the simulations could be performed. Assum- 
ing that the compositions of Table 10 are representative of hydrocarbon 
exhaust emissions for most aircraft engines, the average composition of 
exhaust is 57 percent alkane, 22 percent alkene, 16 percent aromatic, and 5 
percent aldehydes. From the values of Table 10 an average molecular weight 
of 89 g mole"^ can be determined. This molecular weight is close to that of 
a six-carbon molecule, so a six-carbon molecule was used in the computer 
simulation. 
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Figures 2 to 11 show the ozone Isopleths for Mixtures 1 and 3 listed in 
Table 24. Ozone concentration wjrves are given in intervals of 0.1 ppm 
(solid lines) or 0.01 ppm (dashed lines). These isopleths show that Mixture 
1» which contained the highest concentration of alkenes, was the most reac- 
tive producer of ozone. Mixture 3, which contained the lowest percentage of 
alkenes, was the least reactive producer of ozone. Of the different mix- 
tures, Mixture 3 is the mixture most representative of typical aircraft 
mixtures. However, the data base did not seem sufficient to warrant a 
study relating airplane population to ozone production. 

The isopleths in general shw a distinctive evolutionary pattern; the 
region of maximum ozone moves toward the upper , ight corner with time. A 
diagonal line along the ridge of any isopleth divides it into two regions. 
Varying the hydrocarbon (HC) concentration at a constant concentration 
would have a significant effect on ozone concentration in the lower right 
region of an isopleth. Varying HC would have a minimal effect in the upper 
left region. Conversely, varying the NO concentration at a constant hydro- 

X 

carbon concentration would have a more pronounced effect in the upper left 
region than in the lower right region. The kinetic mechanism has never been 
validated at hydrocarbon concentrations greater than 18 ppm (as carbon). 
Reactions that are insignificant at lower hydrocarbon concentrations may 
become important at the higher concentrations. Therefore, there is some 
uncertainty in the region from 18 ppm (as carbon) to 30 ppm (as carbon) 
hydrocarbon concentration. 

Isopleths of one and eight hour simulations of Mixtures 2 and 4 are 
presented in Figures 12 to 15. Since Mixtures 2 and 4 contain approximately 
the same amount of alkenes, their isopleths are very similar. Mixture 2 is 
slightly more reactive since it contained a little more alkene than Mixture 4. 


Hydrocarbon— ppm as carbon 








Tota.1 Hydrocarbon--ppm «s carbon 
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GURE 12. OZONE ISOPLETH FOR A SIMULA 
AFTER A ONE-HOUR PERIO 




Total Hydrocarbon— ppm as carbon 




Total Hydrocarbon — ppm as carbon 








49 


C. COMPUTER SIMULATIONS WITH AUTOMOBILE EMISSIONS 

The reactions of a hydrocarbon mixture typical of automobile emissions 
were simulated for comparison with the results of the aircraft mixture simula- 
tions. The Bureau of Mines (1973) performed an extensive study of hydrocar- 
bon and NO^ emissions from ten 1970-1973 automobiles. The results of their 
study are presented in Table 25. The hydrocarbon distributions are presented 
as weight percent of total hydrocarbon. Converting from weight percent to 
moles of carbon for the 40 percent aromatic fuel, one finds that the C 2 -Cg 
group is most representative of paraffins, the C^-Cg group is most repre- 
sentative of the olefins, and the C^+ group is most representative of the 
aromatics. Therefore, in the computer simulations, butane is used for the 
alkane reactions, propylene for the alkene reactions, and toluene for the 
aromatic reactions. The percentages of the alkane, alkene, and aromatic 
groups from Table 25 are approximately 20 percent, 45 percent, and 35 per- 
cent, respectively. Computer simulations were performed using hydrocarbon 
and initial concentrations in Figure 1 and a NO/NO 2 ratio of 9:1. The 
ozone isopleths obtained from these simulations are shown in Figures 16 and 
17. The hydrocarbon mixture from automotive emissions is very similar to 
Mixture 1 of the aircraft mixtures. A comparison of the ozone isopleths of 
Mixture 1 with the ozone isopleths of the automobile mixture shows similar 
trends in ozone behavior. Also, both of these mixtures show approximately 
the same hydrocarbon reactivity upon comparing the isopleths at any given 
time. 


Table 25 


HYDROCARBON DISTRIBUTION IN AUTOMOBILE EXHAUST 
(In weloht oercent) 


tot^i Hydrocarbons 
Automobile (qrams/mtle) 

W2 Old^mobile 
98 with a 45S* 
CIO fn0ine 

1.79 

1971 Ford CaUxle 
with a 3S1-C10 
Engine 

4.88 

1971 Plyiiouth 
Fury III with 
a 360-CI0 Engine 

3.22 

1872 Ford Torino 
with a 35UCID 
Engine 

2.16 

1970 Chevrolet 
Impala with a 
35O-CI0 Engine 

3. S3 

1970 Pontiac with 
a 400-C10 Engine 

S.86 

1970 Volks wagon 
with a 1.600‘CC 
Engine 

2.44 

1971 Chevrolet Vega 
with a 2*300-CC 
Engine 

4.92 

1973 Ford Torino 
with a 3S1>CID 
Engine 

2.79 

1973 Chevrolet 
Impala with a 
350-ClD Engine 

2,16 


Pdrafftns Olefins 


Methane 

Cj-C5 

V 

Ethylene 


3.0 

8.6 

17.4 

9.6 

9.0 

8.7 

7.4 

12.6 

8.4 

9.6 

3.6 

8.3 

12. S 

10.6 

12,7 

10.9 

7.7 

8.9 

10.8 

8.4 

6.3 

7.0 

11.4 

10.0 

10.4 

3.7 

10.1 

21.7 

8.0 

12,8 

6.6 

7.3 

12.8 

9.4 

8.6 

7.4 

9.4 

19.9 

7.9 

8.8 


8.1 

11, S 

11.0 

13.1 

6.1 

8.8 

14.3 

10.6 

10.3 


Ar ofl^tlcs 


Esl 

Benzene 

bL 

Acetylenes 

2.1 

4.5 

39.4 

6.2 

1.S 

4.8 

35.2 

11.8 

1.2 

5.1 

40.1 

5.9 

0.8 

6.0 

31.2 

15.3 

1.1 

5.5 

40.6 

7.7 

3.1 

3.9 

21.7 

5.0 

1.1 

5,2 

39.8 

9.2 

2.8 

3.8 

30.9 

9.1 

1.2 

5,3 

38.5 

6.6 

1.6 

5.2 

36.7 

6.4 


Notes; Fucl-*typ1c«l II!; ^0 percent arof^atk. Data welcjhtcd In accordance 

with the 197S Federal test procedure. All tests conducted at 75^ F ambient 
ten^rature. 

Source: Bureau of Mines (1973) 









53 


IV SENSITIVITY RUNS AND EVALUATION 
OF THE KINETIC MECHANISM 


As with all theoretical formulations, some uncertainties are inherent 
in the kinetic mechanism. These uncertainties may be due to the inaccuracy 
of the rate constants or the omission of important reactions. Minimization 
of these uncertainties is a major goal of the experimentalist. Yet, the 
evaluation of the kinetic mechanism for the reliability of predictions must 
be left to others. 

The Hecht, Seinfeld, and Dodge general kinetic mechanism was evaluated 
by Hecht et al. (1974a). They found that predictions of the kinetic mechanism 
are strongly dependent on the production of radicals (HO^, 0», OH*, and 
RO*). Without these radicals to initiate the reactions, oxidant levels will 
be very low. The following reactions are the main sources of radicals: 

NO 2 + hv -► NO + 0* 

Aldehyde + hv (1 + 6)H02 + 0 - «)R02 + CO 
HNO 2 + hv -*■ NO + OH* 

Unless a small concentration of aldehydes is used as an initial condition, 
the photolysis of aldehydes cannot occur until after they are formed through 
reactions involving hydrocarbons. In all the computer simulations previously 
mentioned, the initial aldehyde concentrations were assumed to be zero. How- 
ever, aldehydes may compose up to 15 percent of the total hydrocarbons emitted 
from aircraft engines. To investigate the effects of the presence of alde- 
hydes at the start of the simulation, we perfonned computer simulations with 
initial concentrations of aldehydes varying from zero to 0.15 ppm. The 
results of these simulations are shown in Figure 18. The simulations with 
initial aldehyde concentrations show high amounts of ozone during the early 



the ozone 

ppm C 
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hours of the simulations. The simulation with Initial aldehyde concentration 
equal to 0.15 ppm shows a more rapid ozone formation In the early hours than 
the simulation with Initial aldehyde concentration at 0.05 ppm. But by the 
end of the simulations the ozone concentrations for the different Initial 
aldehyde concentrations are essentially the same. Therefore, Increasing the 
Initial aldehyde concentration has an effect on the Induction period, but 
little effect on the concentration of ozone near the end of the simulation 
for the Initial hydrocarbon and NO^ concentrations used In these sinnilatlons. 

As In the case of aldehydes, unless some HNO 2 Is assumed to be present 
Initially, photolysis of HNO 2 cannot occur until It Is formed In the heter- 
ogeneous surface reaction: 

NO + NO 2 + H 2 O 2 HNO 2 

Surface reactions Involving and water are strongly dependent on the 
characteristics of the smog chamber. Durbin et al. (1975) performed sensi- 
tivity runs In which they varied the rate constants of the first two reactions 
listed as surface reactions In Table 23. They also varied the Initial HNO 2 
concentrations. The results of their simulations are presented in Figures 
19 and 20. Increasing the rate constant for the formation of HNO 2 causes a 
rapid buildup of ozone in the first few hours of simulation. Yet, the con- 
centrations for the different runs are the same at the end of the simula- 
★ 

tions. The same effects were observed when the simulations were run with 
Increased initial concentrations. The result of a simulation using the 
generalized kinetic mechanism without the heterogeneous HNO 2 chemistry is 
shown in Figure 21. There is a decrease In ozone concentration in the early 
hours, but the final ozone concentration is the same as in the simulation 
with the heterogeneous HNO 2 chemistry. A second simulation was performed 
with an initial aldehyde concentration equal to 0.05 ppm and without the HNO 2 
chemistry. The results (Figure 21) show an increase in ozone formation in 
the early hours and a final concentration similar to that in the simulation 
with no initial aldehyde concentration, but with the HNO 2 chemistry. 

♦ 

Note that these simulations were carried out for only a six-hour period. 
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The rate of the NO 2 photolysis reaction Is dependent upon the Initial 
concentration of NOg. Due to the wide range of values for the ratio of NO 
to NOg In aircraft emissions, sensitivity runs were performed on the general- 
ized kinetic mechanism with varying Initial NO and NO 2 concentrations. Runs 
were made with NO to NO 2 ratios of 2:1 and 20:1. The results of these simu- 
lations are shown In Figure 22 along with those of a simulation using the 
standard NO/NO 2 ratio of 9:1. These simulations show more rapid Initial ozone 
formation with decreasing NO/NO 2 ratios, but changes In the ratio do not 
affect the ozone behavior near the end of the s.^ulations. Therefore, for 
the range of NO/NO 2 ratios that typify aircraft emissions data, the value of 
the ratio chosen (9 in this case) will affect the ozone behavior only in the 
early hours of the simulation. 

A. SENSITIVITY OF REACTIONS IN THE KINETIC MECHANISM 

Sensitivity analyses for each reaction in the kinetic mechanism have 
been performed by Hecht et al. (1974a). Because only small changes have 
been made in the kinetic mechanism (additions or deletions of certain reac- 
tions and son« revision of rate constants), a brief sensitivity study was 
performed on only the previously determined nine most uncertain reactions in 
the kinetic mechanism The reactions considered in this study are shown in 
Table 26. 

The procedure for the sensitivity study was essentially the same as that 
used by Hecht et al. (1974a). The standard concentration-time profile of 
ozone was obtained by computer simulation of Mixture 3 from Table 24, with 
all rate constants at their standard values (Table 23). Then simulations 
were performed in which the rate constant of one of the nine reactions men- 
tioned above was doubled while all other rate constants were unchanged. The 
concentration-time profiles of ozone obtained from these computer simulations 
are based on the same initial conditions as the "standard" simulation, with 
the exception of a single rate constant change. These profiles were then 
compared with the "standard" profile. The areas of the absolute differences 
between these profiles and the standard were calculated. Then the procedure 
was repeated with each rate constant in turn decreased by a factor of 2. 
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Table 26 

REACTIONS CONSIDERED IN THE SENSITIVITY STUDY 


Reacti on 

NO 2 + hv NO + 0* 

O 3 + NO NOg + O 2 
HO^ + NO + NO 2 + OH* 

OH* + NO 2 HNO 3 

HNO 2 + hv NO + OH* 

Aldehyde + hv •*> Products 
Aldehyde + hv (1+6 )H0^ + (1-6)R0^ + CO 
Alkene + OH* Aldehyde + aHOg + (l-a)RO^ 
NO + NO 2 + H 2 O 2 HNO 2 


Rate Constant 
0.25 min“^ 

25.2 ppm’^ min"^ 

800 ppm“^ min’^ 

4 -1 -1 

1 X 10 ppm min 

1.75 X 10’^ min‘^ 

1.13 X 10"^ min"^ 

1.13 X 10"^ min’^ 

4 -1 - 

2.5 X 10 ppm min 

-7 -2 

2.6 X 10 ppm min 


The areas of absolute differences for each rate constant were averaged. 

These averages were normalized to the area under the curve of the "standard" 
profile and multiplied by 100 to give the sensitivity due to a change by a 
factor of 2 in the rate constants. The results of this sensitivity study are 
shown in Table 27. Note that the two aldehyde reactions have been combined 
into one reaction in the sensitivity study. 

In addition to the sensitivity, the uncertainty of each rate constant 
has an effect on predictions made with the kinetic mechanism. If a rate 
constant is well established, little uncertainty is introduced into the pre- 
dictions by its use, even though the results may be highly sensitive to that 
rate constant. On the other hand, a rate constant with a high degree of un- 
certainty may have little effect on predictions If the mechanism is relatively 
insensitive to that constant. Table 28 lists the uncertainty factors of 
the eight reactions. Hecht et a1. (1974a) postulated an ad hoc index called 
the S*U index to tabulate the major sources of uncertainty in the kinetic 
mechanism. S*U index is defined as: 

S*U = Sensitivity x Uncertainty Factor 
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Table 27 

SENSITIVITY OF THE REACTIONS 


Rank 

1 

2 

3 

4 

5 

6 

7 

8 


Reaction Sensitivity 

O 3 + NO NOg + O 2 34 

NOg + hv NO + 0. 29 

NO + HO 2 NO 2 + OH* 16 

NOg + OH* HNO 3 14 

Aldehyde + h\> Products 11 

Alkene + OH* -»• Aldehyde + 6 

oHO^ + (1 -q)R02 

HNO 2 + hv NO + OH* 3 

NO + NO 2 + HgO 2 HNO 2 2 


Table 28 

UNCERTAINTY FACTORS OF THE REACTIONS 

Uncertainty 



Reaction 

Factor* 

O 3 + NO NO 2 + O 2 

1.3+ 

NO 2 + hv 
NO + HOj 

NO + 0* - 
NO 2 + OH* 

1.4 

3.2^ 

NOg + OH- 

• -V HNO 3 

2 . 0 '^ 

Aldehyde 

+ hv Products 

3 

Alkene + 

OH* Aldehyde + 

1.2 

oHO^ + 

(l-alROj 

HNO 2 + hv -> NO + OH* 

3.0 

NO + NO 2 

+ H 2 O 2 HNO 2 

10.0 


★ 

Subjective estimates, 
t 

Hanipson and Garvin (1975). 


Tiffin 
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Table 29 lists the reactions studied in order of decreasing S*U index. 
The reaction of NO with HO^ radicals shows the highest S*U index, while the 
reaction of alkenes with OH* radicals shows the least S*U index of this 
group. The uncertainties in the predictions of the kinetic mechanism may 
be reduced if more reliable rate constants are found for the reactions 
investigated in this study. 


B. EXTENSION OF THE KINETIC MECHANISM TO INCLUDE LONGER CHAIN HYDROCARBONS 

The kinetic r/iechanism has been validated only for butane and propylene. 
Therefore, extending the mechanism to include longer chain hydrocarbons 
involves, at a minimum, performing sensitivity studies of the parameters 
that are intended to account for the longer chain hydrocarbons. Conkle et al. 


Table 29 

COMBINED SENSITIVITY AND UNCERTAINTY OF THE REACTIONS 


Rank 

Reaction 

S*U Index 

1 

NO + HOg NOj + OH* 

51.2 

2 

O 3 + NO ■*■ NO 2 + O 2 

44.2 

3 

NO 2 + hv -*■ NO + 0 

40.6 

4 

Aldehyde + hv ■*■ Products 

33.0 

5 

NO 2 + OH* -*• HNO 3 

28.0 

6 

NO + N0» + H 2 O -*■ 2 HNO 2 

20.0 

7 

HNO 2 + hv NO + OH* 

9.C 

8 

Alkene + OH* -* Aldehyde + qHO^ + (l-alRO^ 

7.8 
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(1975) observed that n-ocUne, propylene, 1-butene, and p-xylene are the most 
abundant molecules emitted from a T-56 con^ustor. We performed sensitivity 
runs on different combinations of these molecules with aircraft exhaust 
emissions consisting of 55 percent alkanes, 25 percent alkenes, and 20 percent 
aromatics (Mixture 3). Three points on the hydrocarbon/NO block (Figure 1) 

A 

were used as starting conditions. The first two simulations represent situa- 
tions where exhaust having high hydrocarbon and low concentrations is 
emitted into the atmosphere (the hydro carbon/NO^^ ratio is high). The first 
simulation was performed with a hydrocarbon/NO^ ratio of 30:1 (ppm as carbon/ 
ppm as NO^). The starting conditions were 3.0 ppm C of hydrocarbon and 0.1 
ppm of NO . The second simulation was performed with a hydrocarbon/NO 

A A 

ratio of 20:1. The starting concentrations were 6.0 ppm C for the hydro- 
carbons and 0.3 ppm of NO . The third simulation, starting with hydrocarbon 

A 

and NO concentrations equal to 3.0 ppm C and 1.2 ppm, respectively, is 

A 

representative of systems with low hydro'-arbon emissions and high NO^ emis- 
sions. The new rate constants for these simulations are shown in Table 30. 


Table 30 

RATE CONSTANTS FOR LONGER CHAIN HYDROCARBONS 


Rate Constants 

^R^sctjion (p pm~^ min ~^) 

Alkane (n-octane) + 0* ->• RO 2 + OH* 500 

Alkane (n-octane) + OH* ->• RO^ + H 2 O 1 x 10^ 

Aromatic (p-xylene) + 0* RO^ + OH* 290 

Aromatic (p-xylene) + OH* -»■ RO^ + H 2 O 2.7 x 10^ 

Alkene (1 -butene) + OH* Aldehyde + oHO^ + (1 - a)R02 6 x 10^ 


Sources: Doyle et al. (1975), Hampson and Garvin (1975), Project Clean 

Air (1970). 
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Again, the computer simulations were performed under conditions similar 
to those of smog chamber experiments (I.e., constant Irradiation during the 
simulation period and inclusion of surface reactions). The results of the 
three simulations are shown In Figures 23, 24, and 25. For the two simula- 
tions representing systems with high HC/NO ratios (Figures 23 and 24), 

A 

changing to the longer chain hydrocarbons does not have a large effect on the 
ozone concentration compared to the simulation with a lower HC/NO^^ ratio 
(Figure 25). In Figures 23 and 24, the longer chain molecule combinations 
Increase the ozone concentration much faster In the early hours. At any 
given time there Is a greater concentration of ozone than In the shorter 
chain system. The two different longer chain systems give similar ozone con- 
centrations, which is to be expected, because propylene and 1 -butene differ 
by only one carbon atom. Since the longer chain hydrocarbon system Is more 
representative of typical hydrocarbons emitted from aircraft engines, the 
ozone isopleths for longer chain systems as a function of time would show 
a faster buildup of ozone in the early hours for the lower hydrocarbon and 
NO concentrations. In Figure 25, changing to the longer chain hydrocarbons 

A 

has a big effect on the ozone behavior for the low HC/NO^ system. The final 
ozone concentration is more than a factor of two higher in the longer chain 
systems than in the shorter chain system. 

Besides varying the specie^ in the sensitivity runs, we also varied the 
closure parameters, 6 and e. As discussed previously, these closure parameters 
are dependent on the specific alkane and alkene molecules used in the kinetic 
mechanism. For the mixture consisting of n-octane, 1 -butene, and p-xylene, 
the closure parameters were varied from their original values of 0.5 to the 
new values of 0.25. The closure parameters, 6 and e, are functions of the 
length of the hydrocarbon chain. These parameters limit the number of RO^ 
and HO 2 radicals formed in the following reactions: 

Aldehyde + hv -*• (1 + 6 )H 02 + (1 - 6 )R 02 CO 

Aldehyde + OH* -+ 6 HO 2 + (1 - 6 )R 02 

RO 2 + NO >• Aldehyde + 0.5(1 + e)li 02 + 0.5(1 + e)C0 + 0.5(1 - clRO^ + NO 2 












-octane 








I'jtane, propylene. Toluene 

i-octane. propylene, p-xy1ene 

-octane, 1-butene, p-xylene 

-octane, 1-butene, p-xylene (« ■ e - 0.25) 
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Consequently, the closure parameters limit the amount of ozone formed. If 
the amount of formaldehyde Is a substantial fraction of the total aldehydes, 
the first two reactions will be less Important. Also, If a substantial amount 
of ROg Is In the form of pero;o^methyl (CH^Op radicals, then the third reac- 
tion will also be less Important. With longer chain molecules, the fractions 
of formaldehydes and CH^O^ are expected to be lower than for shorter chain 
molecules. Therefore, for the longer chain molecule systems lower values of 
6 and e are expected. The value of 0.25 for s and c may be a lower limit 
because the rate at which the longer hydrocarbons break down to smaller hydro- 
carbons (shorter lengths) as a function of chain length Is presently unknown. 
For exanple. If the rate of breakdown of the chain Is fast, then the amounts 
of CHjOg and formaldehyde may be larger than the amounts determined by the 
fractions given by 6 and e. Therefore, 6 and e are estimated to have a lower 
limit of 0.25 for the longer chain hydrocarbon systems in this study. 

Figures 23, 24, and 25 show the results of the closure parameter changes 
at the three different starting conditions. Again, for the simulations with 
high HC/NO^ ratios (Figures 23 and 24), varying the closure parameters has 
little effect on the ozone concentration when predictions are compared with 
those of the simulation of the system with the low HC/NO ratio (Figure 25). 

A 

Varying the closure parameters for the low HC/NO^ system resulted in twice 
as much ozone at the end of the simulation. 

The overall effect of the presence of longer chain hydrocarbons (includ- 
ing changes in the closure parameters) is an increase in ozone concentration 
of approximately 20 to 30 percent at the end of the simulation period for 
the system.*: having high HC/NO ratios (Figures 23 and 24) and an increase in 

A 

ozone concentration of approximately a factor of 6 at the end of the simula- 
tion period for the shorter chain system for the low HC/NO^^ system (Figure 25), 
The unreactive low hydrocarbon and high NO system was affected more by the 

A 

longer chain hydrocarbons than was the more reactive system (high hydrocarbon 
and low NO^^) . 


iiir Ilf 
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V SIMULATIONS OF MIXED AUTOMOBILE AND AIRCRAFT EMISSIONS 


In the previous chapters, we presented emissions data for jet aircraft 
and for autotnobiles. Each of these sources emits pollutants that can 
lead to ozone production in the presence of ultraviolet radiation. In this 
chapter we discuss simulations of systems containing both jet emissions 
and automotive emissions. Such systems might occur in the real world where 
an airport is located in or near a region of significant automotive emissions, 
such as an urban area or a freeway. We are interested in these systems 
because of the possibility of enhanced ozone production— that is, the possi- 
bility that a system with mixed emissions will produce more ozone than a 
system wherein jet and automotive emissions are kept segregated. 

The creation and validation of a complete numerical model of a mixed 
system at any particular airport would have to include the effects of such 
factors as: 

> Transport and dispersion of pollutants in the atmosphere as a 
function of space and time, 

> The spatial and temporal character of emissions. 

> The dilution of pollutants with clean air. 

> The extent of mixing of jet and automotive emissions. 

> Variations in the duration and intensity of sunlight. 

> The reactions of pollutants before and during mixing. 

> Differences in the types of hydrocarbons emitted from 
different sources. 

> The relative concentrations of pollutants. 

> Induction effects. 

A model of this complexity is beyond the scope of this study. Consequently, 
instead of real physical situations, we discuss several idealized physico- 
chemical models of the reaction and mixing processes, with the intent to 
illustrate rather than document the phenomena that may occur when jet and 
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automotive emissions mix. These models were constructed to meet three 
objectives: ease of calculation, ease of interpretation, and physical 
significance. The effects of dilution and other variations in pollutant 
concentrations were not included in the models. Although Inclusion 
of these effects causes only slight changes in the total amount of 020 ne 
formed from a given quantity of pollutants, it can cause large spatial 
variations in ozone concentrations, which would obscure enfancement effects. 

Before discussing the models and simulations of idealized mixed systems, 
we present a simple example that, although not explicitly related to a 
physical situation, will provide insight into the causes of enhanced 
ozone p'Jduction. 


A. THE POSSIBILITY OP ENHANCED OZONE PRODUCTION FROM MIXING- 
A SIMPLE EXAMPLE USING INITIAL CONDITIONS 

To illustrate the possibility of enhanced ozone production, we 
plotted isopleths of the maximum 1-hour average ozone concentrations 
for Mixture 3 (Figure 26). (Plots of maximum 1-hour average isopleths, 
are more convenient for the present discussion because they do not depend 
on time.) In Figure 26, Line A, at an HC/NO^ ratio of l.’u (ppm as C/ppm 
as NO^), represents typical aircraft emissions, which have high hydro- 
carbon concentrations and low NO^ concentrations. Line B, at an HC/NO^^ 
ratio of 3.0 is representative of automobile emissions. The isopleth plot 
shows that the initial concentrations of Point 1 on Line A (12 ppmC 
hydrocarbons and 0.1 ppm NO ) produce a maximum 1-hour average ozone 

A 

concentration of approximately 0.25 ppm. The initial concentrations of 
Point 2 (3 ppmC hydrocarbons and 1.0 ppm NO^) produce a maximum 1-hour- 
averago ozone concentration of approximately 0.35 ppm. Now, consider the 
averaging of the initial concentrations of Points 1 and 2. This may be 
visualized as the combination and mixing of two air parcels of equal 
volume, each parcel containing one set of initial concentrations of HC. and 
NOj^. The average of the initial concentrations of Points 1 and 2 is 
75 ppmC hydrocarbons and 0.55 ppm NO^, which is shown as Point 3 on 
Line C in Figure 26. One might expect that the maximum 1 -hour-average 
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ozone concentration produced from the average initial concentrations of 
HC and NO^ (Point 3). would equal the average of the ozone concentrations 
at Points 1 and 2, v/hich is 0.3 ppm. But this is not the case; the ozone 
concentration at Point 3 is about 0.47 ppm. The difference of 0.17 ppm 
ozone is the "enhanced" ozone production. 

The cause of the enhanced ozone production may be explained in terms 
of the HC/NO ratios of the emissions. The isopleth diagram shows that for 

X 

a given concentration of HC, there is an HC/NO^ ratio for which ozone 
production is a maximum. The same is true for a given concentration of 
NO , although the HC/NO ratio is different. The HC/NO ratios that maximize 

X X X 

ozone production witli respect to HC and with respect to NO both have 

A 

numerical values between the HC/NO ratio of jet emissions on the one hand 

X 

and automotive emissions on the other. Thus, we expect that enhanced 
ozone production will often occur when jet and automotive emissions mix. 

The simple example described above omits many possibly important 
influences on ozone production such as mixing, dilution, and variations in 
sunlight. It also assumes that the reactivities of the hydrocarbons in 
jet and automotive emissions are equal. As noted above, accounting for these 
influences in a complete numerical model is beyond the scope of this study. 

It is possible, however, to examine ozone enhancement more carefully by 
using models more complex than the simple example just described. The 
following sections describe three such models, and the results of simulations 
using these models. 


B. DESCRIPTION AND PHYSICAL SIGNIFICANCE OF THE MODELS 
1 • Model l--Emissions Mixed Initially 


The first model is that in which jet and automotive emissions mix 
before any chemical reactions occur. This model includes two cases (see 
Figure 27). In the first case, the initial concentrations of pollutants 
are simply added; this represents injection of jet and automotive emissions 


Air Parcel Containing 
Emissions from Source A 



Source B 


Emissions from Source B 


a) Case 1— Addition of emissions from Source B to an air parcel containing 
emissions from Source A. It is assumed that no chemical reaction 
occurs until emissions from both sources are mixed homogeneously 
within the parcel. 



Air Parcel of Volume X 
Containing Emissions 
from Source A 



Containing Emissions 
from Source B 



Air Parcel of Volume 2X 
Containing Emissions 
from Both Sources 


Time = 0 Time > 0 

b) Case 2--Mixing of emissions without reduction in volume, or "averaging" 
of emissions over two parcels. Again, it is assumed that no chemical 
reaction occurs until emissions from both sources are mixed homogeneously 
within the parcel. 


FIGURE 27. 


MODEL 1 , CASES I AND 2— EMISSIONS MIXED INITIALLY 


Into the same air parcel. In the second case, the concentrations of 
pollutants are averaged; this represents the mixing, at some distance from 
the sources of emissions, of two separate air parcels carrying different 
emissions. The two cases of Model 1 approximate physical situations in 
which, during the daytime, mixing is very rapid compared with the rate of 
chemical reaction, such as when air containing jet emissions mixes with 
air from a nearby freeway. Model 1 also represents the physical situation 
during nighttime, when no photochemical reactions are occurring, and air 
parcels containing different types of emissions have a long time (and, 
depending on meteorological conditions, a large volume) in which to mix. 

In terms of the concentrations of pollutants after mixing, most physical 
situations probably lie between the two cases of Model 1. 

2. Model 2--Em1ssions Injected into an Air Parcel Containing 
Reacting Emissions 

Model 1 and the base case (unmixed emissions in separate air parcels) 
represent the extremes of transport limited and reaction limited situations. 

If the production of photochemical smog were a linear pracess, those extremes 
would provide us with all the information necessary to predict ozone/precursor 
relationships for combined automotive and aircraft emission systems. 

However, because of the nonlinearities of photochemical smog production, mixing 
in a partly reacted system may show phenom.ena not seen in fully mixed 
reacting systems. 

In the second idealized model, emissions are injected into an air 
parcel containing other emissions that are already undergoing chemical 
reaction (see Figure 28). As one example, this model represents an air 
parcel passing over an airport and then a freeway, if, during the travel 
time between the airport and the freeway, the jet emissions in the air parcel 
undergo significant chemical reaction. In using Model 2, we considered 
both jet einissions injected into automotive emissions and vice-ve»'sa. 

Although Model 2 is physically similar to Case 1 of Model 1, it differs 
chemically because of the assumption of continuous chemical reaction. 


Air Parcel Containing 
Emissions from Source A That 
Are Undergoing Chemical Reaction 



Source B 


Emissions from Source B 


FIGURE 28. MODEL 2— EMISSIONS INJECTED INTO AN AIR PARCEL 
CONTAINING REACTING EMISSIONS 


3. Model 3--Mixinq of Two Air Parcels Containing Reacting Emissions 


In the third idealized model, an air parcel containing, say, reacting 
jet emissions mixes slov^'ly with an air parcel containing reacting automotive 
emissions (Figure 29). It is assumed that no new emissions are added 
during the mixing. Model 3 is physically similar to Case 2 of Model 1, 
but differs chemically because of the continuous chemical reaction before 
and during mixing. 

In addition to these three models, we also performed simulations with 
varying UV irradiation during the simulations, and with averaged rate 
constants and closure parameters for jet and automotive emissions. The 
effects of these changes are presented in a later section. 

C. INITIAL CONDITIONS FOR SIMULATIONS OF MIXED SYSTEMS 

The computer simulations discussed in the remainder of this chapter 
are based on mixtures of hydrocarbons and oxides of nitrogen chosen to 
represent airport and freeway emissions. This section presents details of 
the selection of the following items: 


Slow Mixing 



Air Parcel of Volume X 
Containing Reacting 
Emissions from Source A 



Air Parcel of Volume X 
Containing Reacting 
Emissions from Source B 


Air Parcel of Volume 2X 
Reacting Emissions from 
Source A 

Reacting Emissions from 

R 


Time = 0 


Time > 0 


FIGURE 29. MODEL 3-MIXING OF AIR PARCELS CONTAINING REACTING EMISSIONS 


> Hydrocarbon/NO ratios representative of airport and 

A 

freeway emissions. 

> Hydrocarbon compositions representative of airport and 
freeway emissions. 

> A ratio of airport emissions to freev/ay emissions. 

> An initial concentration of hydrocarbons appropriate for 
use in the simulations. 

> Variation of intensity of sunlight during the simulations. 

In these simulations we used an HC/NO ratio of 31.0 (ppm as carbon/ 

A 

ppm as NO^) to represent the total airport emissions, in contrast to the 
average ratio of approximately 50 that was presented in Section II-D for 
jet emissions alone. Tne lower ratio takes into account the emissions of 
access and service vehicles, which operate at slower speeds and produce 
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emissions having higher HC/NO ratios than automobiles driven at highway 

A 

speeds (as shown by the EPA data in Table 22). A hydrocaroon/NO^ ratio of 
3.0 (ppm as C/ppm as NO ) was used for the freeway emissions. 

A 

A hydrocarbon composition representative of airport emissions was derived 
primarily from the data of Conkle et al. (1975). The composition used was 
57 percent alkanes, 18 percent alkenes, and 25 percent aromatics. Actually, 
Conkle et al. reported 7 percent aromatics and 17 percent aldehydes. The 
aldehyde value may be high; Lozano et al. (1968) and Groth and Robertson 
(1974) reported low values (Section II-B). Also, Lozano et al. (1968) and 
Groth and Robertson (1974) reported that aromatics were approximately 25 
to 30 percent of the total hydrocarbons. Therefore, a hydrocarbon composi- 
tion consisting of 25 percent aromatics was used for the simulations. Initial 
aldehyde concentrations were not included in view of the insensitivity of 
the ozone behavior to variations in these concentrations after the induction 
period (Section IV-A). We calculated an "averaged" hydrocarbon from data 
reported by Conkle et al. (1975) in order to convert hydrocarbon emissions 
into units of moles; similarly, an "averaged" molecule was used to convert 
ppm NO 2 into moles of NO^. The "averaged" hydrocarbon was a Cg molecule 
(hexane). Therefore, the molecular hydrocarbon/NO ratio used in the simu- 

A 

lations was 3.0 (see Appendix). 

A similar procedure was used to find the hydrocarbon composition for 
freeway emissions. The Bureau of Mines (1973) reported a detailed breakdown 
of hydrocarbon compositions from automotive emissions (Table 25). Based on 
these data, it was determined that a hydrocarbon composition consisting of 
20 percent alkanes, 45 percent alkenes, and 35 percent aromatics would be 
suitable for use. The representative molecules for each hydrocarbon group 
are butane for alkanes, propylene for alkenes, and p-xylene for aromatics. 

The "averaged" hydrocarbon calculated from Table 25 was a C^ molecule. This 
gives a value of approximately 0.4 for the molecular hydrocr,rbon/N0^ ratio 
in units of ppm as an "averaged" molecule per ppm as NO (see Appendix). 

A 
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In most of the simulations we assumed that the concentrations of hydro- 
carbons in air parcels containing airport and freeway emissions are approx- 
imately equal. To ccsnpare this assumption with "real world” conditions, we 
estimated the emissions from San Francisco Airport and from the Bayshore 
Freeway, which is adjacent to that airport. From traffic count data for 
automobiles on the freeway, data on the number of flights to and from the 
airport, and the emissions data reviewed ir Chapter II, we found that the 
hydrocarbons emitted fron the airport equal the hydrocarbons emitted from 
30 miles of the freeway. Because jet emissions have a higher HC/NO ratio 
than automotive emissions, the emissions of nitrogen oxides from the airport 
are only about 10 percent of the emissions from the 30 miles of freeway. 
Airport emissions are generally more concentrated at first than are freeway 
emissions. But an airport is similar to a point source of pollutants, and 
consequently its emissions often undergo more rapid dilution than freeway 
emissions, which come from a line source. It may be expected that at some 
time during the dilution of jet emissions and entrainment of freev/ay emis- 
sions, the conditions in an air parcel approximate our assumed conditions. 

The total initial concentration of hydrocarbons used in each simulation 
was chosen so that the resulting ozone concentration would lie in a range 
between the federal standard of 0.08 ppm and 0.5 ppm, which is typical of 
a very smoggy day in the l.os Angeles Basin. Another reason for our choice 
of initial concentrations was to study airport emissions and freeway emis- 
sions that have similar potential for ozone production. If the ozone con- 
centrations produced separately by the airport emissions and by the freeway 
emissions were widely different, the simulation results would be highly 
influenced by dilution. As noted above, this dilution might obscure enhance- 
ment effects. 

In all the computer simulations presented in the following sections, we 
attempted to simulate physical systems more realistically by using a varying 
NO 2 photolysis rate constant that corresponds to the variation in intensity of 
sunlight on an equinox day near SS'' latitude (Jeffries et al . , 1974). The max- 
imum NO 2 photolysis rats, and in most cases the maximum ozone concentration. 
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Table 31 

INITIAL CONDITIONS FOR SIMULATIONS 


Initial Condition 

Airport Emissions 

Freeway Emissions 

Hydrocarbon Concentration 
(ppm as C) 



Alkanes* 

0.0855 

0.05 

A1 kenes* 

0.027 

0.1125 

Aromatics* 

0.0375 

0.0875 

Total HC 

0.15 

0.25 

NO Concentration (ppm as NO ) 

A s\ 



NO 

0.045 

0.54 

NOg 

0.005 

0.06 

Total 

0.05 

0.60 

HC/NO Ratio (ppm as "averaged" 
molecule/ppin as NO^) 

3.0 

0.42 


♦In simulations in which hydrocarbons were handled explicitly, rather than 
as generalized or "lumped” species, alkenes were represented by propylene, 
aromatics by p-xylene, and alkanes by n~octane in the airport emissions and by 
butane in the freeway emissions. 



occurred at six hours (or noontime) In the computer simulations. The var- 
iation In sunlight* and hence in UV Irradiation, affects the following 
reactions In the chemical kinetic mechanism: 


Reaction Rate Constant 

NOg + hv -*■ NO + 0* 

HgOg + hv 20H* 

Aldehyde + hv Products 
Aldehyde + hv (1 + 6)H0^ + (1 - «)R0^ + CO 
HNOg + hv NO + OH* 

Table 31 lists the Initial conditions for the computer simulations. A 
simulation was performed with each HC/NO^^ ratio to determine the amount of 
oxidant (ozone) produced. Then, a combination of the freeway ratio with 
the aircraft ratio was used in a set of simulations. In all simulations, 
the initial ratio of NO to NOg In both airport and freeway emissions was 
assumed to be 9. This assumed value Is close to measured values for aircraft 
emissions (Table 7). Measurements of ambient NO concentrations show that 

A 

between 6 a.m. and 9 a.m. on weekdays, the NO/NOg ratio In the Los Angeles 
Basin ranges from approximately 2.0 to 10.0 (Roberts and Roth, 1971). Thus 
the ratio of NO to NOg for automotive emissions was chosen to be 9.0. 

D. SIMULATION RESULTS FOR INITIALLY MIXED SYSTEMS 

The Initial conditions developed In the previous section were usea In 
simulations of photochemical smog formation. This section presents the 
results of simulations for systems In which airport and freeway emissions 
are well mixed before chemical reactions begin. Also presented are the 
results of tests performed to determine whether our parameter estimation 
methods Influenced the ozone enhancement effects seen in the simulations. 


0.0036 k^ 
0.0045 k^ 
0.0045 k^ 
0.07 k, 
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In all simulations a modified version of the Hecht, Seinfeld, and 
Dodge (1974) kinetic mechanism was used, and the physical processes of 
transport and dispersion were not taken into account. Instead, we used the 
models presented in Section B. We wish to emphasize these limitations 
because turbulent mixing generally has a significant effect on photochemical 
smog formation in the atmosphere. 

!• Model 1 Case 1— Simulation Results 


As described above, in Model 1 Case 1 airport and freeway emissions 
are added to a single air parcel. Thus the initial hydrocarbon concentra- 
tion is the sum of the initial hydrocarbon concentrations presented in 
Table 31 for airport emissions (0.15 ppm C) and for the freeway emissions 
(0.25 ppm C), or 0.40 ppm C. In the simulation of Model 1 Case 1 the alkanes 
were treated explicitly in the chemical kinetic mechanism as described in 
Chapter IV. For this case, one might expect that the ozone concentration 
in the absence of enhancement would be the average of the ozone concentra- 
tions generated by twice the airport emissions in a single air parcel 
unmixed with freeway emissions and by twice the freeway emissions unmixed 
with airport emissions. But Figure 30 shows that the addition of airport 
and freeway emissions generates much more ozone than either twice the air- 
port emissions taken alone or twice the freeway emissions taken alone. 

The difference between the ozone concentration predicted for Model 1 Case 1 
and the average of the other two ozone concentrations in Figure 30 is the 
enhanced ozone production. 

The ozone concentrations shown in Figure 30 decrease after about 600 
minutes of simulation because of the decreasing NO 2 photolysis rate, which, 
in turn, is caused by ^he decreasing intensity of sunlight. This process 
can be clarified by examination of the following reactions. 
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^ u ■ ■ 

NOg + hv 


NO + 0- 

0 ) 


0 + Og + M 


Oj + M 

( 2 ) 

' ( 

O 3 + NO 


CM 

0 

CM 

0 

ss: 

(3) 


03 + NO 2 


NO 3 ♦ 02 

(4) 


These are the main reactions for formation and destruction of ozone. The 
rate of formation of "odd" oxygen can be written as follows: 


kjCNO^] -yoito^iw 

d[0,] 

-jf- ■ kjCtlJCOjKH] -KjtNOJCOj] -k^CNOjKOj] 

d[0j d([0] + £0,]) , . 

-3T " '■ dt ■ - (¥"03 * k^[N02])[O3] 

Making the steady-state assumption for the oxygen atoms (I.e. , d[0]/dt « 0), 
we obtain 

^ • kiCNOj] - (kjENO] . [O 3 ]) 

Making a second steady-state assumption for ozone (i.e., d[0j]/dt “ 0), 
the ozone concentration as a function of Is found; 

r« 1 

“ kpW^nqiTRjp’ 



I 


! 

i 

I 
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For a situation in which the value of (the NOg photolysis rate constant) 
is kept constant, the ozone concentration can always build up as the con- 
centration of NO decreases. When k^ varies with the sunlight intensity, 
k^ is small in the early hours of the day and the NO concentration is high. 
Therefore, ozone formation is delayed. Near the end of the day k^ goes to 
zero. 


As k^i decreases, the numerator of the steady-state expression for 
ozone (I.e., the rate of NO 2 photolysis) decreases. The ozone in the system 
reacts with the remaining NO, which lowers the rate kjCNO] until it nearly 
equals k^[N02]. The time at which kj[N0] « k^[N023 depends on the ratio 
of NO to ozone in the system. If the amount of NO is small compared to 

A A 

the amount of ozone, the ozone destruction reactions are slow. Normally 
when k^ is large the reaction of ozone with NO is very much faster than 
ozone destruction by N02^k^[N023]. Then the steady state expression is 
approximately 



k^tNOg] 

“igwor 


j 

In this case the amount of ozone produced is independent of the absolute 
amount of NO^; it depends primarily on the ratio of NO 2 to NO. Hence, 
when irradiation intensity is invariant with time, the ratio of NO 2 to NO 
controls the amount of ozone at all times. However, the ratio of total 
NO^ to ozone becomes Important when solar radiation is varied in the simu- 
lations. For automative emissions the NO /ozone ratio is high, therefore 

A 

the ozone concentration decreases as the sunlight decreases. Thus, the 
steady-state relationship is followed more closely in systems with a high 
ratio of NO to ozone. As k, decreases rapidly the low NO (airport) 

XI X 

system "lags'* behind the steady-state value and ozone may even remain at 
appreciable concentrations after the steady-state value goes to zero (after 
the sun sets), provided there is insufficient NOj^ remaining in the system 
to consume it. 


2 . 


>tode1 1 Case 2*"Sfwu1at1on Results 


a. Alkane Reaction Rate Constants Averaged 

In Nodel 1 Case 2, an air parcel of volume X containing airport emissions 
mixes with an air parcel of vo1un» X containing freeway emissions. It Is 
assinned that both emissions are mixed h(xm>geneous 1y throughout a volume of 
2X before any chemical reaction occurs. The Initial simulation of Model 1 
Case 2 was performed using averaged alkane reaction rate constants In the 
chemical kinetic mechanism, rather than handling butane and n-octane reac* 
tions explicitly. The results of this slimilatlon and of simulations of 
airport emissions alone and fre^ay emissions alone are shown In Figure 31. 
Ozone enhancement Is observed; the predicted ozone concentrations for 
>todel 1 Case 2 are significantly higher than the predicted concentrations 
for the unmixed missions. 

b. Alkane Reaction Rate Constants Treated Fxpllcltly 

The simulation results of Model 1 Case 2 presented in Figure 31 were 
oDtalned assuming an averaged alkane reaction rate constant, and assuming 
that all alkanes produce the same radical (ROp In the kinetic mechanism. 

We performed other simulations of Model 1 Case 2 to determine whether the 
ozone enhancement effect shown in Figure 31 was Influenced by these assump- 
tions. The Initial conditions for these simulations are similar to those 
used In the original simulation of Model 1 Case 2, (Figure 31), with the 
exceptions noted: 

(1) Model 1 Case 2 simulation using "explicit" alkanes (butane 
and n-octane) and their rate constants, and assuming that 
these alkanes produce the same radical (ROp. (The results 
of this simulation are shown as Curve 1 in Figure 32.) 

(2) Repeat of simulation (1) with butane and n-octane producing 
different radicals and aldehydes. (Curve 2 In Figure 32.) 
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(3) Repeat of simulation (2) with closure parameters in 
the kinetic mechanism changed to 0.25 {from 0.5) for 
reactions involving airport emissions. (Curve 3 in Figure 32.) 

The original simulation results for Model 1 Case 2 from Figure 31 are repro- 
duced in Figure 32 to facilitate comparison. The similarity of the curves 
in Figure 32 demonstrates that the assumptions used in the initial simula- 
tion have very little effect on the observed enhancement of ozone production. 

3. Model 2— Simulation Results 


Model 2 is an idealization of a common physical situation: An air 
parcel passes over one emissions source* and then the same air parcel passes 
over another emissions source. Model 2 simulations require specification 
of a number of variables, such as which source emits into the air parcel 
first, how much time elapses before the air parcel passes over the second 
emissions source, and how rapidly the second source injects emissions into 
the air parcel. Each simulation in this section represents an air parcel 
containing either airport or freeway emissions at the concentrations listed 
in Table 31, Other emissions are injected into this parcel at a constant 
rate, beginning at the start of the simulation and continuing until the 
Injected emissions total 0.15 ppm C for airport emissions or 0.25 ppm C for 
freeway emissions. All the simulations were performed with explicit alkanes, 
i.e. butane in freeway emissions and n-octane in airport emissions, but the 
products and closure parameters were not separated. Table 32 lists the 
Model 2 simulations performed and the maximum ozone concentrations in the 
simulations. The results of these simulations are shown in Figures 33 
and 34. The results of the Model 1 Case 1 simulation are reproduced in 
both figures for comparison purposes. 



Table 32 

MODEL 2 SIMULATIONS 
Emissions Present 


Simulation 

Initial ly 

(In Concentrations Listed 
in Table 31 ) 

Injected 

Emissions 

Injett|on 

Time 

(hours) 

Maximum Ozone 
Concentration 
IPPjnL ^ 

1 

Freeway 

Airport 

6 

0.23 

2 

Freeway 

Ai rport 

9 

0.20 

3 

Ai rport 

Ai rport 

6 

0.12 

4 

Airport 

Ai rport 

9 

0.11 

5 

Airport 

Freeway 

4 

0.22 

6 

Ai rport 

Freeway 

10 

0.14 

7 

Freeway 

Freeway 

4 

0.19 

8 

Freeway 

Freeway 

10 

0.11 


■gf — ” — 

In each simulation the rate of injection was adjusted so that the total 
amount of emissions shown in Table 31 was injected at a constant rate 
in the number of hours listed here, after which injection ceased. 


r. 











It is difficult to define or demonstrate ozone enhancement in simula- 
tions with continuous injection of emissions. In such simulations, the 
total mass of pollutants increases with time, unlike all of the simulations 
(including the base cases) discussed in the previous sections. In addition, 
the injection of NO inhibits ozone formation. It is still possible to 
observe enhancement, however, by comparing continuous injection simulations 
with one another. For example, in Simulations 1 and 3, airport emissions 
are injected into reacting freeway and airport emissions, respectively 
(Figure 33). The mixing of the two types of emissions in Simulation 1 pro- 
duces a maximum ozone concentration of 0.23 ppm; Simulation 3, with only 
airport emissions, produces 0.12 ppm ozone. It may be objected that this 
difference is not due to enhancement, but rather to the greater mass of 
pollutants in Simulation 1 than in Simulation 3 at any given time. (The 
initial conditions are 0.25 ppmC in Simulation 1 and 0.15 ppmC in Simula- 
tion 3.) This objection cannot be raised agaii t the enhancement observed 
by comparing Simulations 5 and 7, with freeway emissions injected into reactive 
airport and freeway emissions, respectively (Figure 34). Although the amounts 
and rates of injection of freeway emissions are the same for Simulations 5 
and 7, tiie pollutant mass (ppmC) in Simulation 7 always is greater than 
in Simulation 5 because of the difference in initial conditions: Simulation 

7 initially has 0.25 ppmC freeway emissions, and Simulation 5 has 0.15 ppmC 
airport emissions. Despite this difference, the maximum ozone concentration 
in Simulation 5 is 0.22 ppm, compared to 0.19 ppm fnr Simulation 7. A Similar 
enhancement effect may be noted by comparing Simulations 6 and 8. 

The inhibition of ozone formation by NO can also be demonstrated through 
simulations with continuous injection of freeway emissions. (As shown in 
Table 31, it is assumed that freeway emissions contain much more NO than 
airport emissions.) The inhibition by NO may be seen by comparing Simula- 
tions 7 and 8 with the results of freeway emissions alone (Figure 31). In 
Simulations 7 and 8, freeway emissions are being injected into reacting 
freeway emissions over four hours and ten hours, respectively. By the end 


of ^he Injections Simulations 7 and 8 both contain twice as much total pol- 
lutant (In ppmC) as the simulation of freeway emissions alone. But the 
maxlmtan ozone concentrations for Simulations 7 and 8 are 0.19 and 0.11 ppm» 
respectively, compared with 0.19 ppm for freeway emissions alone. Note also 
that, although the maximum ozone concentrations reached In Simulation 7 and 
In the simulation of freeway emissions alone are the same, the rate of . 
ozone formation is different. In the latter simulation ozone begins to 
accumulate after about 200 minutes of simulation time, but In Simulation 7 
ozone does not accumulate until after Injection stops at 240 minutes. In 
Simulation 8, the injection of NO over a ten-hour period Inhibits ozone 
formation through most of the simulation, resulting In the low maximum ozone 
concentration. 

4. Model 3— Simulation Results 


Model 3 Is an Idealization of a physical situation In which two air 
parcels, each containing reacting emissions, slowly merge with one another. 
In the Model 3 simulation, one air parcel contained freeway emissions and 
the other airport emissions, and the two air parcels were assumed to mix at 
a rate of 5 percent per hour. The ozone maximum in each parcel was slightly 
greater than 0.22 ppm and occurred near 600 minutes (Figure 35). At that 
time diffusion was only about 40 percent complete, yet the enhancement phe- 
nomenon Is apparent from comparison with the ozone maxima produced by free- 
way emissions alone and airport emissions alone (0.20 and 0.15 ppm ozone 
respectively, from Figure 31). 

The predicted ozone concentrations in the two parcels of Model 3 are 
generally not as great as those in the initially mixed system of Model 1 
Case 2. There is one exception: For a brief time in the morning, airport 

emissions produce a faster ozone formation rate than does the initially 
mixed system. This faster rate early in the simulation occurs because the 
low concentration of NO in airport emissions is oxidized to NO 2 in a short 
time, so ozone formation begins early. 
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E. CONCLUSIONS 

We have emphasized the nonlinear nature of ozone formation kinetics and 
the requirements for caution In generalizing about such a process. However* 
the Idealized models recounted In this chapter do have some claim to physical 
significance, and in each of them we observe some degree of enhancement of 
ozone fomatlon caused by a mixing of airport and freeway anisslons. Thus, 
our analysis leads us to conclude that such enhanc«nent occurs regularly In 
the atmosphere. The simulations presented in this chapter were formulated 
to demonstrate enhancement effects: enhancement In the atmosphere will be 
insignificant In many cases due to the mixing of air parcels having very dif- 
ferent concentrations of pollutants or other forms of dilution. 


VI COIICLUSIOIIS 


The purpose of this research is to provide a basis for estimating 
the contribution of jet exhaust to smog. We began by searching the 
available literature on jet exhaust emissions in an attempt to determine 
both their composition and quantity. Virtually all the hydrocarbon 
emissions ivere found to come from the taxi-idle mode. Various engines, 
fuels, and operating conditions were found to give a range of hydrocarbon 
mixes. By computer simulation, these mixes were shown to be capable of 
generating high concentrations of ozone In the presence of nitrogen oxides 
and sunlight. The computer simulations used hydrocarbon emissions data 
and various amounts of nitrogen oxides as starting conditions, as in a 
smog chamber experiment. The results of the simulations were qualitatively 
similar: each mixture of hydrocarbons showed two regions of roughly equal 
potential for ozone formation, one at high HC/NO ratios and the other at 

A 

low ratios. The maximum potential for ozone formation occurred at an 
intermediate HC/NO ratio. 

A 

The sensitivity of these simulations to various parameters was studied 
and the relevance of these computer simulations to the actual atroospher-e 
was discussed. In general, a change in an input parameter caused one of 
two effects: either the induction time to the ozone peak changed or the 
maximum amount of ozone formed was affected. The length of the induction 
period was most sensitive to the presence of photochemical radical precursors 
(such as nitrous acid and aldehydes), uncertainties associated with sensitive 
reactions in the kinetic mechanism, and the UO/NO^ ratio. The parameters 
that had the greatest effect on the amount of ozone formed were the closure 
parameters (which have not yet been established using data from smog chanter 
experinents) and the rate constants that are intended to make the kinetic 
mechanism applicable to longer chain hydrocarbons. Surprisingly, the anxsunt 
of ozone formed was sensitive to the values of the closure parameters and rate 
constants only at low HC/NO^^ ratios. Published jet emissions data show that 
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these emissions have high HC/NO ratios. Automobile emissions, on the other 
hand, have low HC/NO^ ratios, but the original closure parameters and rate 
constants were validated for these emissions at low HC/NO^ ratios. 

Our results indicate that atmospheric blending of jet emissions (high 
HC/NO ratio) with automobile emissions (low HC/NO ratio) may produce a 

* A 

mixture having an intermediate HC/NO^^ ratio associated with increased ozone 
production. Control strategies for ozone are usually most effective when 
they shift the HC/NO emissions ratio away fr<»ii this intermediate ratio. 

A 

For the case of jet emissions, however, the possibility of their mixing with 
autonKJbilc exhaust must be considered. The ozone isopleths from our simula- 
tions indicate that reductions in nitrogen oxide emissions alone would pro- 
duce the largest reductions in ozone for the jet emissions by themselves. 

But this strategy would be effective only for airports located in areas with 

it 

no sources of NO^^ ; it would not control the enhancement of ozone formation 
that might arise from the mixing of jet emissions with nitrogen oxides from 
other sources (such as automobiles or fossil-fueled power plants). 

According to Bahr and Gleason (1974), the application of existing tech 
nology could reduce hydrocarbons emitted by jets during the taxi-idle mode 
by a factor of ten. Such a reduction in aircraft hydrocarbon emissions would 
have a number of effects. The overall emissions (and likewise the overall 
potential Tor oxidant formation) would be reduced. The HC/NO ratio of jet 

X 

emissions would more closely match typical urban ratios, which are due mostly 
to automobile emissions. A sufficient change in the HC/NO^ ratio of jet 
emissions would eliminate the possibility of an enhancement effect from the 
mixing of automobile and aircraft emissions. If the HC/NO^ emissions ratio 
for jet aircraft was reduced to that of automobiles, future studies like the 
present one would I'equire experiiiHsntal validation of th:"* closure parameters 
for long chain molecules as previously mentioned. The iiydrocarbon mixtures 
used In the simulations illicit also require revision, since reducing the 
HC/NO ratio of an engine's emissions might change tiie composition of its 

A 


Rasmussen (197i?) has shown that rural areas have higli HC/NO^ ratios. 
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exhaust hydrocarbons. And lastly, such a reduction might reduce the rela- 
tive importance of the taxi- idle mode as a source of emissions. 

It is important to keep in mind that the conclusions drawn from this 
study are based on results derived for a v/ell mixed system; chemistry was 
examined in depth but transport and mixing were not. If future planning 
needs warrant more careful and accurate exan.Jnation of the impact of jet 
emissions, it will be necessary from an analytic standpoint to incorporate 
the kinetic mechanism in a mesoscale air pollution simulation model that 
describes the relevant atmospheric processes. Results derived from exer- 
cising a model of this type would lead to less restrictive, more general 
conclusions. 


Appendix 

CONVERSION FACTORS 

Due to the different methods of reporting emissions data (e.g., lbs., 
lbs. /hour, lbs,/1000 lbs. fuel, ppm as carbon, etc.), a more uniform system- 
of I'its was used in this report: All data taken from studies by different 
investigators were reported in the units they used, but all hydrocarbon/NO^ 
ratios were reported as moles of carbon/moles of NO 2 . All data of hydro- 
carbon emissions were converted to moles of carbon from pounds by multiplying 
by 37.88 moles as carbon/lb. All data of oxides of nitrogen emissions (un- 
less specified otherwise) were converted from units of lbs. to moles of 
NO 2 by multiplying by 9.85 moles as N02/lb. 

The sections on synerg.stic effects required emissions per total LTD 
cycle. Therefore, where emission rates were reported for individual modes 
(in units of lbs. /hour or kg/hour), the rates were multiplied by the fraction 
of time spent in each mode. The time factors are taken from Table 3. 


The input data for the computer simulations require concentrations 
in units of ppm as "some" molecule. Therefore, a converting factor is needed 
to convert ppm as carbon to ppm as a molecule. For aircraft systems, the 
study by Conkle et al . (1975) of the T-56 combustor was the most detailed. 
Even though the T-56 engines constitute only 11.2 percent of the major 
engines in the U.S. Air Force (Table 32), they are used in 15.8 percent of 
the flying hours. Conkle et al. reported the following compositions for 
hydrocarbons emitted at 15 psig combustor pressure: 57 percent alkanes, 

18 percent alkenes, 7 percent aromatics, and 17 percent aldehydes. Multi- 
plying these percentages by the molecular v/eight reported by Conkle et al. 
for each group, an average molecular weight of 86 g/mole is found. The 
hydrocarbon ^6^14 has the same molecular weight. Therefore, this molecule 
is used to convert from noles of carbon to moles as an "average" nwlecule. 

The converting factor is 0.14 moles as CgH.j^/moles as carbon. For the 
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Table 32 

USAF AIRCRAFT ENGINE USAGE^ 


Percentage 
of Major 


Percentage 
of Flying 


Engine 

Engines 

Hours 

J-57 

30.1 

26.3 

TF-33 

9.3 

17.8 

T-56 

11.2 

15.8 

J-85 

10.5 

10.5 

J-79 

15.6 

9.5 

J-69 

6.6 

7.2 

J-60 

1.4 

2.5 

T-76 

.9 

1.6 

TF-30 

2.8 

1.5 

J-33 

2.2 

1.5 

TF-39 

.9 

1.1 

J-75 

2.0 

1.0 

T-58 

1.4 

1.0 

TF-41 

1.7 

.8 

J-65 

1.1 

.5 

T-64 

.5 

.4 

T-S3 

.6 

.3 

T-400 

.5 

.3 

J-71 

.4 

.3 

J-47 

.3 

.1 


Rased on data from AR.C/wrArR for 19,036 installed active 
otKjines for the first (luarter of 1972. 

Source?: Naugle (197'1). 
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automobile system, the Bureau of Mines (1973) reported that most of the 
alkanes emitted were C^-C^ molecules, most of the alkenes were C^-Cg 
molecules, and the aromatics were Cy+ molecules- They also reported that 
hydrocarbon emissions were 20 percent alkanes, 45 percent alkenes, and 35 
percent aromatics. This gave an average molecular weight around 68 g/mole. 
A molecule was assumed to be the average molecule for the automobile 
system. The converting factor is 0.18 moles as C^/moles as carbon. 

NO^ concentrations were calculated by multiplying by 1.44 moles as 
NO^/mole as The average molecular weight of NO^ was found by assuming 

the ratio of NO/NO^ to be equal to 9. 
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